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Using a generalization of the Madelung transformation, we derive the hydrodynamic representa¬ 
tion of the Klein-Gordon-Einstein equations in the weak field limit. We consider a complex self¬ 
interacting scalar field with a A|v?| 4 potential. We study the evolution of the spatially homogeneous 
background in the fluid representation and derive the linearized equations describing the evolution 
of small perturbations in a static and in an expanding universe. We compare the results with 
simplified models in which the gravitational potential is introduced by hand in the Klein-Gordon 
equation, and assumed to satisfy a (generalized) Poisson equation. Nonrelativistic hydrodynamic 
equations based on the Schrodinger-Poisson equations or on the Gross-Pitaevskii-Poisson equations 
are recovered in the limit c —> +oo. We study the evolution of the perturbations in the matter 
era using the nonrelativistic limit of our formalism. Perturbations whose wavelength is below the 
Jeans length oscillate in time while perturbations whose wavelength is above the Jeans length grow 
linearly with the scale factor as in the cold dark matter model. The growth of perturbations in the 
scalar field model is substantially faster than in the cold dark matter model. When the wavelength 
of the perturbations approaches the cosmological horizon (Hubble length), a relativistic treatment 
is mandatory. In that case, we End that relativistic effects attenuate or even prevent the growth of 
perturbations. This paper exposes the general formalism and provides illustrations in simple cases. 

Other applications of our formalism will be considered in companion papers. 
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I. INTRODUCTION 

Scalar fields (SF) play an important role in particle 
physics, astrophysics, and cosmology m- Their evolu¬ 
tion is usually described by the Klein-Gordon (KG) equa¬ 
tion [MU which can be viewed as a relativistic extension 
of the Schrodinger equation 00 . The KG equation was 
one of the first attempts to unify the ideas of quantum 
mechanics and Einstein’s theory of special relativity. 1 It 
describes the evolution of a SF (p(x, t) whose excitations 
are bosonic particles of spin zero. These particles are neu¬ 
tral for a real SF and charged for a complex SF. The KG 
equation usually involves a potential V (ip) that takes self¬ 
interaction into account. For example, 7r-mesons, Higgs 
bosons and axions are described by the KG equation. SF 
are also present in the theory of superstrings as well as in 
a large number of Kaluza-Klein and supergravity models. 

The coupling between the KG equation and gravity 
through the Einstein equations, leading to the Klein- 
Gordon-Einstein (KGE) equations, was first considered 
in the context of boson stars l9liTf2[. 2 Initially, the study 
of boson stars was motivated by the axion field, a pseudo- 
Nambu-Goldstone boson of the Peccei-Quinn phase tran¬ 
sition, that was proposed as a possible solution to the 


1 The KG equation was actually discovered by Schrodinger before 
he found the equation that now bears his name [2]. 

2 Boson stars are described by complex SFs. Self-gravitating sys¬ 
tems described by real SFs are not static but, instead, are pe¬ 
riodic with both the spacetime geometry and the matter field 
oscillating in time. For that reason they are called oscillatons 
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strong CP problem in QCD. In the early works of Kaup 
[9] and Ruffini and Bonazzola [10], it was assumed that 
the bosons have no self-interaction. For a boson mass 
to ~ 1 GeV/c 2 , the maximum mass of boson stars is very 
small, of the order of M max ~ 10~ 19 Mq. This corre¬ 
sponds to mini boson stars like axion black holes. The 
mass of these mini boson stars may be too small to be as- 
trophysically relevant. They could play a role, however, 
if they exist in the universe in abundance or if the axion 
mass is extraordinarily small (less than to ~ 10” lo eV/c 2 ) 
leading to macroscopic objects with a mass M max com¬ 
parable to the mass of the sun (or even larger) 09 ■ it 
has also been proposed that stable boson stars with a 
boson mass to ~ 10 -17 eV/c 2 could mimic supermassive 
black holes (M ~ 10 6 Mq, R ~ 10 7 km) that reside at 
the center of galaxies ESKH3- On the other hand, Colpi 
et al. [14] considered the case where the bosons have 
a repulsive self-interaction described by a A|y>| 4 poten¬ 
tial and showed that, for to ~ lGeV/c 2 , the maximum 
mass of boson stars can be of the order of the solar mass 
Mq, similar to the mass of neutron stars. Therefore, a 
self-interaction can significantly change the physical di¬ 
mensions of boson stars, making them much more astro- 
pliysically interesting. 

It has also been proposed that dark matter (DM) halos 
may be made of a SF described by the KGE equations 
(see, e.g., for recent reviews). Actually, at the 

galactic scale, the Newtonian limit is valid so DM halos 
can be described by the Schrodinger-Poisson (SP) equa¬ 
tions or by the Gross-Pitaevskii-Poisson (GPP) equa¬ 
tions. In that case, the wavefunction if>(x, t) describes 
a Bose-Einstein condensate (BEC) at T = 0, and the 
self-interaction of the bosons is measured by their scat- 
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tering length a s . Therefore, DM halos could be gigantic 
quantum objects made of BECs. The wave properties 
of bosonic DM may stabilize the system against gravi¬ 
tational collapse, providing halo cores and sharply sup¬ 
pressing small-scale linear power. This may solve the 
problems of the cold dark matter (CDM) model such 
as the cusp problem and the missing satellite problem. 
The scalar field dark matter (SFDM) model and the 
BEC dark matter (BECDM) model, also called TDM 
models, have received much attention in the last years. 
In the non-interacting case [38MT] . the mass of the 
bosons must be extremely small, of the order of to ~ 
2.57 x 10 _2O eV/c 2 , in order to account for the mass and 
size of dwarf DM halos that are completely condensed. 
Ultralight scalar fields like axions may have such small 
masses. This corresponds to “fuzzy cold dark matter” 
[45) . On the other hand, when a repulsive self-interaction 
is taken into account [571 ESI EMS], the characteristics 
of dwarf DM halos can be reproduced with a much larger 
boson mass m = 1.69 x 10 _2 eV/c 2 * (satisfying the limit 
in < 1.87 eV/c 2 obtained from cosmological considera¬ 
tions EH) and a scattering length a s = 1.73 x 10 5 fm 
(satisfying the constraint 47ra 2 /m < 1.25 cm 2 /g set by 
the Bullet Cluster m- Dwarf halos are purely con¬ 
densed. Large dark matter halos have a core-halo struc¬ 
ture resulting from gravitational cooling J21] . They are 
made of a solitonic core (BEC) surrounded by a halo of 
scalar radiation. 

Scalar fields have also been introduced in cosmology. 
The phase of inflation in the early universe is usually 
described by some hypothetical SF, called the inflaton, 
with its origin in the quantum fluctuations of the vac¬ 
uum [7511521 . A variety of inflationary models that in¬ 
clude SFs have been proposed (53], and SF are expected 
to play an important role in determining the dynamics of 
the early universe. Inflation is generally considered to be 
a reasonable solution to many of the fundamental prob¬ 
lems within the standard cosmological model. Inspired 
by the analogy with the inflation, some authors have rep¬ 
resented the dark energy responsible for the present ac¬ 
celeration of the Universe by a SF called quintessence 
[841196] . This treatment seems to fit current observations, 
and unlike the cosmological constant A, the SF evolves 
dynamically, leaving distinctive imprints in the Cosmic 
Microwave Background (CMB) and matter power spec¬ 
trum. 

Since DM may be a SF, it is of considerable interest to 
study the cosmological implications of this scenario. The 
cosmological evolution of a spatially homogeneous non¬ 
interacting real SF described by the KGE equations com¬ 
peting with baryonic matter, radiation and dark energy, 
was considered by Matos et al. [53]. They found that real 
SFs display fast oscillations but that, on the mean, they 
reproduce the cosmological predictions of the standard 
ACDM model. The study of perturbations was consid¬ 
ered by Suarez and Matos [97] for a self-interacting real 
SF described by the Klein-Gordon-Poisson (KGP) equa¬ 
tions and by Magana et al. [55; for a non-interacting 


real SF described by the KGE equations. These stud¬ 
ies show that the perturbations can grow in the linear 
regime, leading, in the nonlinear regime, to the forma¬ 
tion of structures corresponding to DM halos. This is 
in agreement with the early work of Khlopov et al. [991 
who studied the Jeans instability of a relativistic SF in a 
static background. 

The case of a complex self-interacting SF representing 
BECDM was considered by Chavanis {57[ 1001 in the con¬ 
text of Newtonian cosmology. His study is based on the 
GPP equations. The Jeans instability of a homogeneous 
self-gravitating BEC in a static background is studied 
in m ■ The evolution of perturbations of BECDM in 
an expanding Einstein-deSitter (EdS) universe is consid¬ 
ered in [1001 . It is found that the perturbations grow 
faster in BECDM as compared to ACDM. Harko m 
and Chavanis 1100] independently developed a relativistic 
BEC cosmology by assuming that the equation of state 
of BECDM is given by P = 27t a s h 2 e 2 /m?c A 7 which cor¬ 
responds to the classical equation of state of BECs [ 102] 
where p is replaced by e/c 2 . However, the identification 
of the energy density e with the rest-mass density p is 
only valid in a weakly relativistic regime, so the extrapo¬ 
lation of their results to early times is not correct (see the 
discussion in [ 103) ). Recently, Li et al. [104] developed an 
exact relativistic cosmology for a complex self-interacting 
SF/BEC based on the KGE equations. They studied the 
evolution of the homogeneous background and showed 
that the Universe undergoes three successive phases: a 
stiff matter era, followed by a radiation era due to the 
SF (that exists only for self-interacting SFs), and finally 
a matter era similar to CDM. They compared their theo¬ 
retical results with observations in order to constrain the 
parameters (to, A) of the SF. 

Instead of working directly in terms of a SF, we can 
adopt a fluid approach and work with hydrodynamic 
equations. In the case of the Schrodinger equation, this 
hydrodynamic approach was introduced by Madelung 
[105] . He showed that the Schrodinger equation is equiv¬ 
alent to the Euler equations for an irrotational fluid with 
an additional quantum potential 3 arising from the finite 
value of h. This hydrodynamic formulation was criti¬ 
cized, or disregarded, by many authors (e.g., Pauli) in the 
early years of quantum mechanics because it lacks a clear 
physical interpretation in the case where the Schrodinger 
equation describes just one particle. However, it takes 
more sense when the Schrodinger equation (or the GP 
equation) describes a BEC made of many particles in 
the same quantum state [102] . In that case, the BEC 
can be interpreted as a real fluid described by quantum 


3 The results of Madelung Esa were rediscovered by Bohm ES] SO 

the quantum potential is sometimes called the Bohm potential. 

The formulation of classical and relativistic quantum mechan¬ 
ics in terms of hydrodynamic equations was also developed by 
Takabayasi 110711108| . 
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Euler equations. 4 This hydrodynamic representation has 
been used by Bohmer and Harko m and by Chavanis 
[5711551 1100) among others in the case of BECDM and in 
the case of BEC stars [35]. This hydrodynamic approach 
has been generalized by Suarez and Matos m in the 
context of the KG equation. They used it to study the 
formation of structures in the Universe, assuming that 
DM is in the form of a fundamental SF with a A (p 4 po¬ 
tential E2- They also studied the phase transition of a 
real SF due to a Zi symmetry of its potential ra- 

In the works mm, the SF is taken to be real and 
the gravitational potential is introduced by hand in the 
KG equation, and assumed to be determined by the clas¬ 
sical Poisson equation where the source is the rest-mass 
density p. This leads to the KGP equations. However, 
this treatment is not self-consistent since it combines rel¬ 
ativistic and nonrelativistic equations. In this article, 
we derive the hydrodynamic representation of a complex 
SF coupled to gravity through the Einstein equations in 
the weak field approximation. In this way, we develop 
a self-consistent relativistic treatment. Throughout this 
work, we use the Newtonian gauge which takes into ac¬ 
count metric perturbations up to first order. We consider 
only scalar perturbations. This is sufficient if we are in¬ 
terested in calculating observational consequences of the 
SF dynamics in the linear regime. After having derived 
the hydrodynamic equations, we study the evolution of 
the homogeneous background and the evolution of small 
perturbations in a static and in an expanding universe us¬ 
ing the hydrodynamic representation. We compare our 
results with those obtained from the heuristic KGP equa¬ 
tions. We argue that this simplified model is not suffi¬ 
cient to study the evolution of the perturbations in the 
linear relativistic regime. 

The main purpose of the present paper is to develop a 
general formalism. Applications of this formalism will be 
considered in following papers. However, for illustration, 
we already present a few applications of our formalism in 
simple cases. 

The paper is organized as follows. In Sec. [TTJ we 
present the theoretical background and introduce the 
main equations of the the paper that are the KG equa¬ 
tion and the Einstein field equations. This allows us to 
set the notations. In Sec. |III[ we consider the KGE 
equations in the weak field limit. We first write the KG 
equation and the Einstein equations in the Newtonian 
gauge. Then, we transform the KGE equations into the 
Gross-Pitaevskii-Einstein (GPE) equations to facilitate 
the connection with the nonrelativistic limit c —> +oo. 
Finally, we write the GPE equations in the form of hy¬ 
drodynamic equations. In Sec. m we consider the case 
of a static background and study the evolution of small 


4 One interesting aspects of BECs is related to their superfluid 
properties. Their velocity field is irrotational but there may exist 
vortical motion due to singular point vortices with quantized 
circulation h/m 11091 [1101 . 


perturbations in an infinite homogeneous Universe (Jeans 
problem). We determine the energy of the homogeneous 
SF, the dispersion relation of the perturbations, and the 
Jeans length. In Sec. [Yj we consider the evolution of the 
Universe induced by a homogeneous SF. We recover from 
the fluid equations the three phases previously obtained 
by Li et al. m-- a stiff matter phase, a radiation phase, 
and a matter phase. In Sec. |VI[ we consider the evolu¬ 
tion of the perturbations in an expanding homogeneous 
Universe. We provide the exact set of linearized equa¬ 
tions, then propose a closed approximate equation for 
the density contrast that has the correct Jeans length. 
In Sec. |VII| we study the evolution of the perturbations 
in the matter era. We show that the nonrelativistic limit 
can be used when the wavelength of the perturbations 
is much smaller than the cosmological horizon (Hubble 
length). Perturbations whose wavelength is below the 
Jeans length oscillate in time while perturbations whose 
wavelength is above the Jeans length grow linearly with 
the scale factor as in the CDM model. The growth of 
perturbations in the scalar field model is substantially 
faster than in the CDM model. When the wavelength of 
the perturbations approaches the cosmological horizon, 
general relativity attenuates or even prevents the growth 
of perturbations. The Appendices present useful com¬ 
plements. In Appendix [A] we determine the constant 
needed to transform the KG equation into the GP equa¬ 
tion. In Appendix[B] we derive the hydrostatic equations 
describing spatially inhomogeneous relativistic SF/BEC 
clusters. In Appendix 0 we recover from our formal¬ 
ism the nonrelativistic equations obtained in the frame¬ 
work of the GPP equations. In Appendix [D] we intro¬ 
duce the generalized KGP equations in which the gravita¬ 
tional potential is introduced by hand in the KG equation 
in a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) 
background space-time, and assumed to be given by a 
generalized Poisson equation where the source is the en¬ 
ergy density e. We show that this model corresponds to 
the limit <b/c 2 —> +00 of the KGE equations. In Appen¬ 
dices [EjG] we regroup all the elements necessary to make 
the numerical applications needed in the paper. A sum¬ 
mary of our results can be found in the Proceedings of 
the X Mexican School on Gravitation and Mathematical 
Physics ra¬ 


il. THEORETICAL BACKGROUND 

A. The Lagrangian of the scalar field 

We assume that DM can be described by a complex 
SF 5 which is a continuous function of space and time 


5 There are several reasons for considering a complex rather than 
a real SF EE|. Complex SFs have been invoked in many dif¬ 
ferent sectors of elementary particle physics in relation to the 
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defined at each point by p{x fl ) = p(x, y, z , t). The action 
of the relativistic SF is 


S v = 


1 


d 4 X\f zz gC lp , 


( 1 ) 


where £ v = £ v (p, p* ,dpp, dp,p*) is the Lagrangian den¬ 
sity and g = det (g^) is the determinant of the metric 
tensor. We adopt the following generic Lagrangian den¬ 
sity 

= M 2 ), (2) 

which is written for a metric signature 
Specifically, we consider a SF potential of the form 

2 2 2 

t r / I iO\ TTb C . . O 777 / . i A , \ 

v(M 2 ) = w M 2 + ^am‘, ( 3) 

where the quadratic term is the rest-mass term and the 
quartic term is a self-interaction term. In view of future 
applications, we shall assume that the SF describes a 
BEC at T = 0 in which all the particles are in the same 
ground state. In that case, the self-interacting constant 
A can be expressed in terms of the scattering length of 
the bosons a s and of their mass m by A = 47 t a s h 2 /m. 
The potential of the SF can be rewritten as 


v{\p\ 2 ) = 


2 h 2 


+ 


2-7T a,m. 


K 2 


(4) 


The self-interaction is repulsive when a s > 0 and attrac¬ 
tive when a s < 0 11021 . The case of a general potential 
of the form V(\p\ 2 ) is treated in | I1 121 . 


B. The Klein-Gordon equation 


The equation of motion for the SF can be obtained 
from the principle of least action. Imposing 6S V = 0 
for arbitrary variations Sip and Sp* , we obtain the Euler- 
Lagrange equation 


D„ 


d£ u 


d{d fl p)* 


d £, 
dp* 


= 0 , 


(5) 


where D is the covariant derivative. For the generic La¬ 
grangian ([2]), this leads to the KG equation 


Op + 2K(M 2 ),^ = 0, 


where □ is the d’Alembertian operator 

□ = D^ v d v ) = -^=d^gg^d v ) 


( 6 ) 


(7) 


Higgs mechanism responsible for mass generation. On the other 
hand, the U(l) symmetry implies the dark matter particle num¬ 
ber (charge) conservation [2]. Finally, complex scalar fields can 
form singular vortices leading to a rich dynamics of the halos 
when they rotate Ell- 


and 


V(\p\ 2 ) 




dV 

d\p\ 2ip ' 


( 8 ) 


For the specific SF potential |3]), the KG equation takes 
the form 


Up + 


h 2 


-p- 


8tt a s m. 

h 2 


\p\ 2 p = o. 


(9) 


The d’Alembertian operator can be written in terms of 
the covariant derivative and of the Christoffel symbols as 


np = D^ v d v p) = g^D^p) 

= g^d il d v p-g^Vl v d a p. ( 10 ) 


C. The energy-momentum tensor 

Taking the variation of the SF action with respect 
to g ^ v , we get 

SS v = l -j d i xV=gT lu ,5g* 1 ', (11) 


where 


2 d£p_ £ 

-L fiis — ^ o uv yi±v / ~-'<p 


dgV v 


( 12 ) 


is the energy-momentum tensor of the SF. For the generic 
Lagrangian ([2]), it takes the form 


Tfiv = -(<9 p,p*d u p + d v p*dp,p) 


9^ 


-g p(T d p p*d a p-V{\p\ 2 ) 


(13) 


By analogy with the energy-momentum tensor of a per¬ 
fect fluid, the energy density and the pressure tensor of 
the SF are defined by e = T® and P/ = — T/ . The con¬ 
servation of the energy-momentum tensor, which results 
from the Noether theorem, writes = 0. 


D. The Einstein equations 

The Einstein-Hilbert action in general relativity is de¬ 
fined by 


Sg ~ 16 nG . 


d 4 a 


■V-gR, 


(14) 


where R is the Ricci scalar and G is Newton’s gravita¬ 
tional constant. Its variation with respect to is given 
by [mj: 


ss g = - 


d 4 xy/-g (r^ -\gpvR^ Sg^, (15) 


167tG 
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where R^ v is the Ricci tensor. The total action (SF + 
gravity) is 5 = S v + S g . The held equations can be 
obtained from the principle of least action. Imposing 
5S = 0 for arbitrary variations in , and using Eqs. 
(11) and (15), we get the Einstein equations 


R- 87rG T 


(16) 


These are a set of 10 equations that describe the funda¬ 
mental interaction between gravity and matter as a result 
of the curvature of space-time. The energy-momentum 
tensor T^ v is the source of the gravitational Held in the 
Einstein Held equations of general relativity in the sense 
that it determines the metric g^ v . The conservation of 
the energy-momentum tensor is automatically included 
in the Einstein equations. 


III. THE KLEIN-GORDON-EINSTEIN 
EQUATIONS IN THE WEAK FIELD LIMIT 


In this paper, we study the KGE equations in the weak 
field limit $/c 2 <C 1. The equations that we derive are 
valid at the order 0($/c 2 ). For <f>/c 2 —> 0, we obtain the 
generalized KGP equations (see Appendix|D|. Of course, 
the limit <f>/c 2 —► 0 is different from the nonrelativistic 
limit c -A +oo leading to the GPP equations (see Ap¬ 
pendix [C]). 


A. The conformal Newtonian Gauge 


We work with the conformal Newtonian gauge Ell¬ 
in general relativity, this gauge is a perturbed form of 
the FLRW line element. We assume that the Universe 
is flat in agreement with the observations of the cosmic 
microwave background (CMB). The general perturbed 
FLRW metric in the comoving frame has the form 


ds 2 = ( 1 + 2— ) c 2 dt 2 — 2a(t)w.icdtdx l 




<t> 

1 - 2^1 5 , 


H, 


dx l dx\ (17) 


where the perturbed quantities ’P/c 2 , 4>/c 2 , Wi, and H ig 
are all <C 1. In this metric, $ represents the gravitational 
potential of classical Newtonian gravity while \I/ is the 
lapse function |I114| . 

We consider the simplest form of the Newtonian gauge, 
only taking into account scalar perturbations which are 
the ones that contribute to the formation of structures 
in cosmology. Vector (which are supposed to be always 
small) vanish during cosmic inflation and tensor contri¬ 
butions (which account for gravitational waves) are ne¬ 
glected |115j . We also assume '!' = $, supposing a uni¬ 
verse without anisotropic stress. Therefore, our line ele¬ 


ment is given by 

1 — 2— j 5ijdx l dx J , 
(18) 

where we recall that <f>/c 2 <C 1. Working with this metric 
enables us to obtain exact equations at the order 0(4>/c 2 ) 
taking into account both relativistic and gravitational 
contributions inside an expanding Universe without hav¬ 
ing to introduce the gravitational potential $ by hand. 



B. The Klein-Gordon equation 


Computing the d’Alembertian ([7]) with the Newtonian 
gauge (fl8]) , we obtain the KG equation 


1 d 2 p 
c 2 dt 2 


3 H dp 
c 2 dt 


1 



A ip 


2 





4 9$ dp 
c 4 dt dt 


(19) 


where H = a/a is the Hubble constant and A = V 2 is 
the usual Laplacian (V is the usual nabla operator). For 
the specific SF potential Q, the KG equation takes the 
form 


1 d 2 p 3Hdp 1 / 4<P\ A 

?~dW + ~di ~ ^ + 1 2 ) * 

4 dp f 24>\ m 2 c 2 

c 4 dt dt \ c 2 ) h 2 ^ 

+ |1 + ^)— 1 ^ = °- 


( 20 ) 


Using the expression (13) of the energy-momentum ten¬ 


sor, the energy density and the pressure are given by 

! / ° - 2 


fl-^ 

dp 

l c 2 J 

dt 


C T °° 2c 2 

+ i( 1+ f) l ^' 2+ nM ! ). 


( 21 ) 




dp 


dt 


-^2 ( 1 +^)|V+| 2 -U(|+| 2 ). (22) 


C. The Einstein equations 

The time-time component of the Einstein equations is 

k -\r= ( 23 ) 
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With the Newtonian conformal gauge, the left hand side 
of Eq. (23) is given by 


r°-Ir = 


3 H 2 




Therefore, the time-time component of the Einstein equa¬ 
tions can be rewritten as 


A<f> 

47rGa 2 


e 

W 


3 H 2 


+ 


3 H 


8t rG 4 ttGc 2 


<9$ 

~dt 


H<& 


(25) 


Using the time-time component of the energy-momentum 
tensor, which represents the energy density e given by Eq. 
(21), we get 


+ 


A<f> 
47rGa 2 
1 


1 

2? 


2a 2 c 2 


2 $ 

1 ~ TT 


dp 


dt 


2<f> \ 2 rrr 


1+ * 


27ra s m. l4 3 H 2 
-\p\ - 


h 2 c 2 


+ 


3H 


8ttG 4ttGc 2 




(26) 


Eqs. (201 and (261 form the KGE equations. This sys¬ 


tem of equations is closed because we consider a universe 
without anisotropic stress (T = <f>). Otherwise, we need 
to write the other components of the Einstein equations. 


which is one of the Friedmann equations [113j . The other 
Friedmann equation is obtained from the Einstein equa¬ 
tion (p5|) that reduces, when = 0, to 


H 2 = 


8t rG 
"3c 2 


£b- 


(32) 


This relation shows that the term —3H 2 /8nG = — e&/c 2 
in the Einstein equation (25) with *!>(x, t) ^ 0 plays the 
role of a neutralizing background. From Eqs. (31) and 
(32), we easily obtain 


a 

a 


4ttG 

lie 2 " 


(eb + 3 Pb)- 


(33) 


E. The Gross-Pitaevskii-Einstein equations 


The KG equation without self-interaction can be 
viewed as a relativistic generalization of the Schrodinger 
equation. Similarly, the KG equation with a self¬ 
interaction can be viewed as a relativistic generalization 
of the GP equation. In order to recover the Schrodinger 
and GP equations in the nonrelativistic limit c —> +oo, 
we make the transformation [2]: 


D. Spatially homogeneous scalar field 


For a spatially homogeneous SF with pb{x,t) = pb(t) 
and t) = 0, the KG equation (19) reduces to 


1 d 2 p b 
c 2 dt 2 


3 H dpb 
c 2 dt 


+ 2V (\p b \ 2 ) ,<p£ = 0. 


(27) 


For the specific SF potential (Q, we explicitly have 


1 d 2 ipb 3 H dpb 
c 2 dt 2 c 2 dt 


2 2 
me 


h 2 


-pb ■ 


&Tta s m 

h 2 


\Pb\ 2 Pb = 0 . 


(28) 


For a homogeneous SF, the energy-momentum tensor is 
diagonal and isotropic, Tjp = diag(et, — Pb, — Pb, — Pb)- 
The energy density e b (t) and the pressure Pb(t) are given 
by 


1 

eb ~ 2d 2 


dpb 

dt 


+ m^i 2 ), 


(29) 


Pb = 


2 c 2 


d<fb 

dt 


-V(\pb\ 2 ). 


(30) 


From these equations, we obtain the continuity equation 


deb 

dt 


+ 3 H{e b + P b ) = 0 


(31) 


p{x,t) = -e~ imc2t ^(x,t). 
m 


(34) 


The prefactor H/m is justified in Appendix |A| Mathe¬ 
matically, we can always make this change of variables. 
However, we emphasize that it is only in the nonrela¬ 
tivistic limit c —► +oo that ip has the interpretation of a 
wave function, and that \ip\ 2 = p has the interpretation 
of a rest-mass density. In the relativistic regime, ip and 
p = \ip\ 2 do not have a clear physical interpretation. We 
will call them “pseudo wave function” and “pseudo rest- 
mass density”. Nevertheless, it is perfectly legitimate to 
work with these variables and, as we shall see, the equa¬ 
tions written in terms of these quantities take relatively 
simple forms that generalize naturally the nonrelativistic 
ones. 


Substituting Eq. (34) in Eqs. (20) and (26), we obtain 


d 0 h 2 d 2 ip 3 S 2 dip 

dt 2me 2 dt 2 2 me 2 dt 

h 2 f 4<f>\ 

+ w ( 1 + ^ 


47t a s h 2 
m 2 



\ip\ 2 ip H —ihHip 


2 h 2 < 9 $ (dip 

me 4 dt \ dt 


ime 2 \ 

—*) 


= 0 , 


(35) 
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A$ 

AttGcl 2 


= (i-l)w 


+ 


h 2 

2 m 2 c i 


/ 2*\ 

dip 

2 h 2 

V c 2 ) 

dt 

2 a 2 m 2 c 2 


1 + 


2 $ 


IW-f 


27r a,fi 2 


M 4 - 


mc z 
3 H 2 




3 H 


8t rG 4 ttGc 2 


r,* 

a +ff,f 


Eq. (351 can be interpreted as a generalized Schrodinger 


equation (in the absence of self-interaction a s = 0) or 
as a generalized GP equation (in the presence of self¬ 
interaction a s / 0). It is coupled to the Einstein equation 

©. 

The energy density and the pressure can be written as 


a pseudo action. Following Madelung, we also define a 
pseudo velocity field as 


v(x, t) = 


VS 


(40) 


where the scale factor a has been introduced in order to 
take into account the expansion of the Universe [97]. We 
note that this velocity field is irrotational. 


(36) 


Substituting Eqs. (39 )-(M0| in the GPE equations (35)- 


(|36|), and separating real and imaginary parts, we obtain 

dS ' 


the system of hydrodynamic equations 


dp TT 1 Id 

- + mp+ - v . (pv) = —- 


"at 


3Hp dS 4p d$ / 2 _ dS 
me 2 dt + me 4 dt V dt 


4$ - 

-— 2 V-{pv), 

ac z 


(41) 


h 2 


2 m 2 c 2 


2 $ 
1 - 2 


dip 


dt 




m \ c 
h 2 


2 a 2 m 2 


dt 


i+5) iwf 


+ r 2 f 1 ^?) w 2 +vm 2 ), 


p = 


h 2 


2 m 2 c 2 


2 $ 
1 ~~ c 2 " 


dip 


dt 




m \ a 
h 2 


6a 2 m 2 


dt 


i+5'i iwf 


jc 2 fi-f) \ip\ 2 -vm 2 ). 


(37) 


(38) 


Eqs. (35) and 


form the GPE equations. In the 
nonrelativistic limit c -7 +oo, we recover the GPP equa¬ 
tions (Cl)-(C2) of Appendix [Cj 


dS[ (VS) 2 
dt 2 ma 2 


K 


■2 ^Vp 

dt 2 


44> 

c 


2 me 2 y/p 

V 24 ^ 


^ ~ 21 1 2ma 2 ^fp mc 2 a 2 


—m<$> — 


47 Ta s h 2 p 


+ J^( d 4\ -13 


2 me 2 \ dt 


„ 2 $ 
z , H 2 

m z \ c z 

4 <9$\ h 2 dp 


c 2 dt ) 4 me 2 p dt ’ 


(42) 


dv 1 ^ h 2 


o 2 +p' 

dt 2 

Vp 


h 2 


2 m 2 a 3 

87t agh 2 ^. 

—- 

3h 2 


4$ 

1 4- 2 " 

c z 


Ay/P 


Vp \ 


-V$-Vp 

a pa 


ac z 


r V(4>u 2 ) + 


-HV — 

4am 2 c 2 \P dt 


1 

2am 2 c 2 
h 2 


dS 

~dt 


v ( 


am*c 


2 ^ 4 v ^ dt p dt 


(43) 


F. The hydrodynamic representation 

Important characteristics of the system are revealed 
by reformulating the KGE equations in the form of hy¬ 
drodynamic equations. This can be done at the level 
of the GPE equations ((351)- ([36]) via the Madelung trans¬ 


formation m ■ To that purpose, we write the pseudo 
wavefunction ip as 


ip(x,t) = ^p(¥Jy s ^/\ 


(39) 


where p = ipip* = \ip\ 2 plays the role of a pseudo rest- 
mass density and S = (1 / 2)ih\n(ip* / ip) plays the role of 


+ 


h 2 


2m 2 c 4 
h 2 


2$ 

1 _ +2" 


2 $ 

~^2a 2 m 2 c 2 l, c 2 
2ira^h 2 


A$ / $ 

47rGa 2 \ c 2 

1 (dp\ 2 P_(dS\ 

4p \dt J h 2 \dt) 


^(Vp) 2 + |(VS) 2 


+ 


1 


rrfic 2 ^ me 2 1 c 2 ) ^ dt 


2$\ dS 


3 H 2 


3 H 


8ttG 4ttGc 2 


d4> 

m +H * 


( 44 ) 






































































where p is a pseudo pressure given by the polytropic 
(quadratic) equation of state 6 


27rctcS 2 


P = 


(45) 


We note that this equation of state coincides with the 
equation of state of a nonrelativistic BEC with a self¬ 
interaction |102j . This coincidence is not obvious because 
Eqs. (41)-(44) are valid in the relativistic regime. The 


interpretation of this equation of state is, however, not 
direct because p and p are a pseudo density and a pseudo 
pressure that coincide with the real density and the real 
pressure of a BEC only in the nonrelativistic limit c —► 
+ 00 . 

The energy density and the pressure can be written in 
terms of hydrodynamic variables as 


h 2 


2 m 2 c 2 

n 2 

2 a 2 m 2 
24- 


24- 
1 - 2 “ 


1 (?P\ ,P_(dS 

4 p \dt J h 2 \ dt 


+- 


1 + 


24- 


1 


.4 p 


(Vp ) 2 


K 2 


(VS ) 5 


- 1 - 


p dS 
m dt 


P = 


h 2 


2 m 2 c 2 


1 - 


24- 


1 

' 2 
1 

4 P 


1 - 


dp 

dt 


24- 


pc 2 + V(p),( 46) 


P_ 

h 2 


dS\‘ 
Ik ) 


h 2 


6 a 2 m 2 


„ 24- 

1 4 —2 


t-(v P f + £(VS)’ 




We note that, in general, the pressure P defined by Eq. 
(47) is different from the pressure p defined by Eq. (451. 


However, we shall find that they coincide for a spatially 
homogeneous SF. 


The hydrodynamic equations (41)-(44) have a clear 
physical interpretation. Eq. (41), corresponding to the 


imaginary part of the GPE equations, is the continuity 
equation. We note that J pd 3 x is not conserved in the 
relativistic regime. However, we will see in Sec. [V] that 
Eq. (41) is consistent with the conservation of the charge 
of a spatially homogeneous SF. Eq. (42), corresponding 


to the real part of the GPE equations, is the Bernoulli or 
Hamilton-Jacobi equation. Eq. (43), obtained by taking 


the gradient of Eq. (42), is the momentum equation. Eq. 
([44]) is the Einstein equation. We stress that the hydro- 


dynamic equations (41)-(44) are equivalent to the GPE 


equations (35)-(36) whic 
the KGE equations ( [20| 

limit c -* + 00 , we recover the quantum Euler-Poisson 


ir are themselves equivalent to 
and (l26|. In the nonrelativistic 


equations (C3)-(C6) of Appendix [C| 


IV. THE CASE OF A STATIC UNIVERSE: 
JEANS-TYPE INSTABILITY 

We first consider the case of a static Universe. This 
amounts to taking H = 0 and a = 1 in the previous 
equations. 


A. The fluid equations 


In a static Universe, the hydrodynamic equations (41 )- 
(44) reduce to 


dp . 1 d ( dS 

a +, ' H = rf#'a 


4 P as ( 2 as 


as (VS ) 2 

dt 2 to 


44- - 

-^V.(pu), 


(48) 


h 2 


d\(P 

dt 2 


2 me 2 sj~p 


+ I U— n~ 


v \r P _ 2 * ( ^ s? 


44- 

■ 2 1 2m y/p 

47t a s H 2 p f 24- 

—mi> - 7 ;— (l-l—— 


2mc 2 


dS_ 

~dt 


m 

2 


h 2 dp 94> 
mc 4 p dt dt ’ 


(49) 


h 2 


dv -> 

s + («.-v)„ = - 2mV 


d 2 Vf>' 
dt 2 

Vp 


2 m 2 


8na s h 2 2 9 . 

- 0 - 5 -V(p4-) - -o V(4>u 2 ) + 


1 + 4J^ 

Vp J 

1 


— V4--Vp 

P 


2 m 2 c 2 
h 2 


v 


dS_ 

~dt 


+ 


h 2 


2 m 2 c 4 


1 -V 


24- 

2 


v ( 

dt p dt 


A4- f 4-' 

4t rG _ ~ c? 1 P 

2 


(50) 


1 

4 p \dt) h 2 \ dt 


h 2 


2$ 
1J —2 

2'ira< < H 2 


l(Vp) 2 + ^(VS) 2 


1 


1 - 


2$\ dS 
Ik' 


The equation of state p(p) associated with a general potential of 
interaction of the form V(|(/?| 2 ) is given in m - 


( 51 ) 
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B. Infinite homogeneous background 

We consider a homogeneous SF at rest extending in 
an infinite space. We have p(x,t) = pb, v b {x,t) = 0, 
and <f 'b(x,t) = 0. We must be careful, however, that the 
phase can depend on time: Sb(x,t) = Sb{t). Actually, 
the equation of continuity ( |48| implies that dSb/dt is a 
constant that we write —E as it represents the energy of 
the homogeneous SF (with the opposite sign). Indeed, 
for a stationary state, the pseudo wavefunction of the 
SF writes ipb(x,t) = yfpbe~ lEt ^ h , so that S b {t) = —Et. 
The energy of the SF is given by the Hamilton-Jacobi 
equation d49b which reduces to 


E 2 


2 me 2 


4na s h 2 p b 

+ E -„- = U, 


yielding 


E = me 


-l+\ 1 + 


8na s H 2 pb 


(52) 


(53) 


H = 0, a = 1, and $6 = 0 by looking for a stationary 
solution of the form ipb{x,t) = y/pbe~ lEt ^ h . Then, using 
Eq. (34}, implying <p b {x,t) = £^ e -i( mc 2 +B)t/R, we 
find that the total energy of the SF, including its rest 
mass, is E tot = E + me 2 . Using Eq. (|53|, we obtain 


E t ot = me a /1 + 


8-7T a s h 2 pb 


(61) 


C. Linear wave equations 

We now slightly perturb the homogeneous SF and con¬ 
sider the evolution of the perturbations in the linear 
regime. This is the relativistic SF generalization of the 
classical Jeans problem m- As usual, we make the 
Jeans swindle m- We write p(x,t) = p b + Sp(x,t ) 
and S(x,t) = —Et + SS(x,t) with Sp(x,t ) -C 1 and 
dSS/dt(x,t) < 1. We also recall that \v(x,t)\ <C 1 and 
<J>(x, t) < 1. Introducing the notations 


pressure of the SF are 


From Eqs. (46) and (47), the energy density and the 

2'Ka,h 2 


£ b — Pb c ~ + 


E 2 


Pb ( 
to \2me 2 


+ E\ + 


-Pb . 


(54) 


Pb 


a = 


1 dSS 

Ik’ 


(62) 


we can write the linearized equations for the perturba¬ 
tions as 


Pb 


= Pb ( 


TO V 2 TOC 2 


E 2 \ 2na s h 2 
+ E — 


2 

3 Pb' 


(55) 


Combining Eq. (|52|) with Eqs. (|54|) and (|55|), we obtain 

! 2 
3 Pb 


, 67 xa s h 2 9 

£b = PbC + 


„ 2’KCLeh 2 9 

P b = - °^pl 56 

m 6 


We note that the pressure P b of a spatially homogeneous 
SF coincides with the pseudo pressure p given by Eq. 
(451. We also note the identities 


1 + 


E 


mc- 


cM 

dt 


1 + 


E 


V • v = 


1 da 
c 2 dt 


4 d$/' 1 E 


dt 


me* 


h 2 d 2 S h 2 


4m 2 c 2 dt 2 


-cjS -1 + 2 


4to 2 


Ad 




(63) 


(64) 


E 


me* 


( E 


1 = 1 


E 


where 


\2 me 2 


c? = P'iPb) = 


2c| 

El ' 


9 ’ 
r .z. 


4-k a s h 2 pb 


(57) 


(58) 


(59) 


is the square of the speed of sound in the homogeneous 
background. In terms of the speed of sound, the energy 
of the SF can be written as 


E 


= -1 


1 + 


2 c 2 


(60) 


Remark: The expression (53) of the energy E of the SF 


can be directly obtained from the GP equation (35) with 


1 + 


E 


me* 


dv 

dt 


n 2 ^ 

4to 2 c 2 I dt 2 


+ 


A<f> 

4-rrGpb 

E 


-c 2 V6- 1 + 2^ V$, (65) 


E 


= (1 + ^U-I (1 


1 + 


C‘ 

— ) 

2 TOC 2 J 


$ 

~2 
C 2 

2^)' (66) 
Taking the time derivative of Eq. (63) and using • 


r 2 E 4> . 

6 - ^^3 I 1 + 


v = A a, we get 
1 d 2 a . 

~di 2 a ~ 1 


E \ d 2 S 4 d 2 4> ( 1 E 


dt 2 


c 2 dt 2 


(67) 
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Taking the Laplacian of Eq. (64), we obtain 

h 2 


1 + 


E 


me* 


A .. A ( K . 1 d 2 S 

Aa ~ 4TO 2 ( c 2 dt 2 

- (1+^)A$. 


— c:AS 


( 68 ) 


Eqs. (67) and ( 68 ), together with Eq. ( 66 ), form a system 


of three coupled wave equations for the evolution of the 
perturbations (6, a , <I>) in the linear regime. 


D. Dispersion relation and Jeans length 


Decomposing the perturbations in Fourier modes 
of the form S(x,t) = 6k exp i(k • r — cot), a(x,t) = 
Ufc exp i(k ■ f— u>t), and <!?(:?, t) = 4+ exp*(fc • r — uit), we 
obtain after straightforward but lenghty calculations the 
dispersion relation 


h 2 


4 m 2 c 4 


w 4 — 


1 + 7 h 2 k 2 3c 2 " 1 

+ 1 + ^ 


3y + 1 2 m 2 c 2 


1 


l + 3 7 


f} 2 l A 

(l-7)^~'yk 2 c 2 + (1^3 1 )k 2 c 


C 
2 2 


bJ 


= o, 


(69) 


where we have introduced the abbreviation 


7 = 


4irGp b 
k 2 c 2 


1 2 c 2 

1 + ^f 

c 2 


(70) 


In the nonrelativistic limit c —► +oo, the dispersion rela¬ 
tion reduces to ES3: 




2 


h 2 k 4 
4 m 2 


4n Gp b . 


(71) 



FIG. 1: Relativistic (solid line) and nonrelativistic (dashed 
line) dispersion relation. Here ko = (16nGm 2 pb/h 2 ) 1 ' 4 . We 
have taken a s = 0 and x = {4nGh 2 pb/m 2 c 4 ) 1 ^ 2 = 0.05. 


Taking u> = 0 in Eq. (69), we find that the Jeans length 


kj is determined by the equation 


4m 2 


c?- 


—4'KGpb ( 1 + -G 2 - 


nh 2 Gpb 


2c 


k 2 j 


1 2 C s 

1 + -T 
C 2 


(72) 


This is a second degree equation in k 2 . Only the solution 
with the sign + has k 2 > 0. In the nonrelativistic limit 


+oo, Eq. (72) reduces to £7) : 


h 2 k 4 

^+f+c 2 fc 2 -4^ = 0 . 


(73) 


Eq. (69) is a second degree equation in w 2 . We can 


show that the discriminant is positive so the relativistic 
dispersion relation has two branches w±(fc). For k -+ 0, 


-(k)^ 


4 to 2 


h 2 


3c 2 
1 + “T 
c 2 


+ k 2 c 2 + 


(74) 


w 2 (£;) ~ —\>k 2 (?. (75) 

O 

For k -+ +oo, 

tJ±(k) ~ k 2 c 2 . (76) 

The branch u\{k) is always positive, corresponding to 
stable (oscillating) modes. The branch w 2 (k) starts from 
0 , decreases, reaches a minimum, increases, vanishes at 
k = kj, and finally tends to +oo (see Fig. [I]). The 
modes are stable (oscillating) for k > kj and unstable 
(growing) for k < kj. However, in the relativistic case, 
the system is stabilized at very large scales since w 2 (k) —> 
0 for k —> 0 while in the nonrelativistic case w 2 ( 0 ) = 
— 4nGpb < 0. This shows that the limits c -+ +oo and 
k —> 0 do not commute. We shall encounter this large- 
scale stabilization again in Secs. [YlB] and [yllGl 


V. COSMOLOGICAL EVOLUTION OF A 
SPATIALLY HOMOGENEOUS SCALAR FIELD 


We consider the evolution of a universe induced solely 
by a spatially homogeneous SF. In the comoving frame, 
we have p(x,t) = p b (t), Vb(x,t) = 0, = 0, and 

Sb{x,t) = Sb{t). We introduce the notation E(t) = 
— dSb/dt which can be considered as the time-dependent 
energy of the spatially homogeneous SF in the comov¬ 
ing frame. The wavefunction o f th e SF is 1/7 (T. t) = 
\J /Ob(t)e —E (*) dt . UsingEq. (34), we have (x, t) = 
m \/ Pb{t)e-W {mc 2 t+J E(t) to) g0 total energy of the 

SF, including its rest mass, is E tot (t) = E(t ) + me 2 . 
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A. Hydrodynamic equations for a homogeneous 
scalar field 


For a spatially homogeneous SF, the hydrodynamic 
equations (41)-(44) reduce to 


dpb 

dt 


ore 1 d ! cn 3-H> 6 

3 Hp b = -j — ( p b E) - -E, 

me at me z 


E 

2 me 2 


+ 1-5 = 


E 2 d ATra s h 2 p b 


2me 2 m 2 

^ 3HH 2 dpb 
4mc 2 pb dt ’ 


3ff 2 _ ft 2 

87tG 2to 2 c 4 


1 / d/Qb 

4p b \ dt 

2ira s h 2 

F Q *-> 




2 E 

Pb 4-yPb 

me- 


(77) 


(78) 


—(79) 


In terms of the total energy E to t(f) = E{t) + me 2 , the 
equation of continuity (771 becomes 


1 dpb 
Pb dt 


1 d-Etot 


3 da 
a dt ' E tot dt 


= 0. 


(80) 


neous SF. 7 Working directly on the homogeneous KG 
equation, Li et al. [104| have shown that a universe filled 
with a relativistic complex SF first undergoes an intrinsic 
stiff matter era, followed by a radiation era due to its self¬ 
interaction, before finally entering in the matter era. The 
stiff matter era occurs when the SF oscillations are slower 
than the Hubble expansion while the radiation and mat¬ 
ter eras occur when the SF oscillations are faster than 
the Hubble expansion. Since we propose in this paper 
to use hydrodynamic equations instead of the KG equa¬ 
tion, it is important to show, for self-consistency, that 
these different regimes can be directly obtained from the 
hydrodynamic equations (78), (79) and (82). One ad¬ 


vantage of the hydrodynamic representation is that we 
do not have to make averagings over the oscillations of 
the SF as in ill) 1 . 


B. Stiff matter era 

At ear ly tim es (a —> 0, t —> 0), we can approximate 
Eqs. (1781) and (l79l by 8 


t2 d dPb 
n dt 2 


3 Hh 2 dpb 


m 2 c 4 yfpb 2m 2 c 4 pb 


dt 


= 0 . 


(85) 


It can be rewritten as a conservation law: 

4( E to tP b a 3 ) = 0 • (81) 

dt 

Therefore, the total energy is exactly given by 

Etot _ Qm . , 

2 3’ 

mcr pba* 

where Q is a constant. This conservation law was found 
by Gu and Hwang [118j directly from the KG equation 
(see also Appendix B of [104j ). It can be shown that 
Q = J J° d 3 x represents the conserved charge density of 
the complex SF. 

The energy density and the pressure of a homogeneous 
SF are 


e b 


h 2 1 

8m 2 c 2 pb 




m 



EpbC 2 + 


2ira s h 2 

m 3 


Pb, 


(83) 


h 2 1 / dp b \ 2 
8 m 2 c 2 pb \ dt ) 



m 


( 1+ 2^) 

2ira s h 2 2 


(84) 


Equations (78), (79) and (82) determine the complete 
evolution of a universe induced by a spatially homoge- 


3 H 2 


h 2 


dpb 


( 86 ) 


8irG 8m 2 c i pb \ dt 
Recalling that H = a/a, Eq. p6] ) can be integrated into 
3m 2 c 4 


Pb 


ATtGh 2 


(— In a) 2 . 


Substituting Eq. (87) in Eq. (85), we obtain 


dH 

dt 


+ 3H 2 = 0 


which yields 


H = —, 
3t 


t 1 / 3 . 


(87) 


( 88 ) 


(89) 


7 A remark may be in order. A SF tpb(t) oscillates in time because 

it has a phase Sb(t). The kinetic and potential energies and 
V(Ph) of a real SF also present oscillations because they keep 
track of the phase Sb(t) of the SF |5311971 [98| . By contrast, the 
kinetic and potential energies \pb\ 2 an d (|| 2 ) of a complex 
SF do not oscillate in time. It is not clear how one can measure 
the oscillations of a complex SF because we do not have directly 
access to field variables such as Pb{t) but rather to hydrodynamic 
variables such as pb(t ), and Pb(t ). In this sense, the hy¬ 

drodynamic representation of the SF may be more physical than 
the KG equation itself. 

8 We can check a posteriori that the terms that have been ne¬ 
glected are indeed negligible when a —>■ 0 and t —>■ 0. 
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Substituting Eq. (89) in Eq. (871 we find that the pseudo 
rest-mass density behaves as 

2 4 

me 

Pb ~ l2irGh 2 


(-In t) 2 . 


(90) 


According to Eq. (82) the energy of the homogeneous SF 
behaves as 


E AirQGh 2 

2 ~ 


1 


1 


oc 


me 


3 me 4 a 3 (— In a ) 2 (—In t) 2 t 


(91) 


Finally, the energy density and the pressure of the ho¬ 
mogeneous SF behave as 


2„2 


Pb = Q) = 


3 H 2 c 


8t tG 24?rGf 2 


1 

oc 


(92) 


These results can be straightforwardly obtained from 
the KG equation in the regime where the SF oscillations 
are slower than the Hubble expansion 11041. Indeed, ne¬ 


glecting the potential energy in Eqs. (29) and (30), we 


get P = e which is the equation of st ate of a stiff fluid 
for which the speed of sound c s = y/P'(e)c is equal to 
the speed of light (c s = c). Substituting this equation 
of state in the Friedmann equations (31) and (32), we 
directly obtain Eq. 


(92). 


Therefore, when t —» 0, a 
relativistic SF behaves as a stiff fluid 1104] , This stiff 
matter era is independent on the scattering length a s of 
the bosons. Therefore, it exists for both self-interacting 
and non-interacting SFs. 


Equations (95) and (96) determine the evolution of the 


Universe induced by the homogeneous SF in the regime 
where its oscillations are faster than the Hubble expan¬ 
sion. 

The energy density and the pressure are given by 


£& — PbC 4- 


6ira s h 2 




Pb, 


2ira s ti 2 2 
Pb 3 Pb • 


(97) 


We note that the pressure Pb{t) of a spatially homoge¬ 
neous scalar field coincides with the pseudo pressure p(t) 
given by Eq. (45). From Eq. (971 we get 


Pb = 


m 3 c 4 
727ra<A 2 


/ 2Aira s h 2 
1 4 _ , £b - 1 


(98) 


This is the equation of state obtained by Colpi et al. m 
in the context of boson stars (see also [32, 46, 104] ). 

We can now determine the asymptotic behavior of the 
homogeneous SF in the regime where its oscillations are 
faster than the Hubble expansion, i.e. in the regime that 
follows the stiff matter era. 

(i) At late times t —► +oo (pb —>• 0, a —► +oo), we get 


Pb 


Qm 

„3 ’ 


E 4 - 7 T a s h 2 Q 


(99) 


3H 2 

8nG 


Pb 


a ~ ( 67 rGQmt 2 ) 1 / 3 , (100) 

a 3 


C. Radiation and matter eras 

At later times, f ollo wing the stiff matter era, we can 
approximate Eqs. (78) and (79) by 9 


3 H 2 
87 tG 


(E +1 \ E= i ^pL, 

\ 2mc z J m z 


2Tra s h 2 2 E f E 

= Pb~\ - 3~2~Pb 4-- ' - 

m J c 


me 2 \ 2 me 2 


(93) 


+ 1 Pb. (94) 


Substituting Eq. (93) in Eq. (94), we obtain 


3 H 2 
8nG 


— Pb 4" 


677 a eh 2 


3 9 Pb ’ 
m^c z 


(95) 


Eq. (931 is the same as Eq. (52) in the static case, so we 


still have the identity (53), where now E and pb depend 


on time. Combining Eqs. (53) and (82) we get 


Pb 14- 


%'Ka s K 2 pb\ 1 ^ Qm 


= T (96) 


9 These equations correspond to the terms that have been ne¬ 
glected in Sec. IV Bland that now become important. 


6b ~ pbC , 


_ 2na tt h 2 0 / ^ x 

p b ~ —sV^b - 0. (101) 


The SF/BEC universe behaves as a pressureless matter 
fluid. 

(ii) At early times t —> 0 (pt, —> 4 - 00 , o —> 0), we get 


Pb 


Q 2 m 5 c 2 \ 1 ^ 3 1 E l'8Tra s h 2 Q\ L/A 1 

87t a s h 2 ) a 2 ’ me 2 \ m 2 c 2 ) a’ 


1/3 


( 102 ) 


3 H 2 6ira s h 2 

8nG m 3 c 2 


2 3 f Q 4 Trma s h 2 \ 1/3 1 

Pb ~ 2 [ -?- ) (103) 


l G 3 Q 4 


sti 2 


1/12 


t 1 / 2 , 


£b 


67 ia 8 h 2 


_ 2 

3 Pb' 


Pb 


3 £b ' 


(104) 


(105) 


The SF/BEC universe behaves like radiation. As em¬ 
phasized by Li et al. [104j . the radiation era is due to 
the self-interaction of the SF (a s ^ 0). There is no such 
phase for a non-interacting SF (a s = 0). 
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Remark: the case of a SF with an attractive self¬ 
interaction (a s < 0) is theoretically interesting because 
it leads to a model of universe without initial Big Bang 
singularity. At t = 0, the universe has a finite scale fac¬ 
tor a* = (12v / 37rQ|a s |S 2 /m 2 c 2 ) 1 / 3 and a finite density 
/?* = m 3 c 2 / 12ir\a s \h 2 . Two evolutions are possible: (i) 
a normal evolution in which the density decreases as the 
scale factor increases until the system enters in the mat¬ 
ter era, and (ii) a phantom-like evolution in which the 
density increases as the scale factor increases until the 
universe enters in a de Sitter era with a constant den¬ 
sity Pmax = m 3 c 2 /8Tr\a s \h 2 . The choice of the evolution 
between these two behaviors cannot be decided a priori. 
This interesting dynamical system will be studied in a 
specific paper (in preparation). 


A<F 


AnGpba 2 


= 1 


E ( E 


me 2 \2mc 2 


+ 1 


s -4 


E 
me 2 
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h 2 


8m 2 c A pl 

$ 


( 

me 2 \2me 2 

9 95 , 


+ 1 )S 


dt 

H 2 <F 1 


c 2 4m 2 c 4 c 2 p 2 


dfb 

dt 


dpb 
dt 
2 


2$ 


3 H 


AirGpbC 2 \ dt 


d4> 


+ H4> 


where 6 , a and c 2 are defined as in Eqs. 
where Pb(t) is now time-dependent. 


(109) 


(59) and (62) 


VI. EVOLUTION OF THE PERTURBATIONS 
IN THE EXPANDING UNIVERSE 


B. Approximate equation for the density contrast 


We now slightly perturb the homogeneous SF and con¬ 
sider the evolution of the perturbations in the linear 
regime. 


A. Linearized equations 

The linearized hydrodynamic equations obtained from 
Eqs. (41)-(44) describing the evolution of small inhomo¬ 


geneities in the expanding Universe are 


1 + 


E 


mc~ 


cM 

dt 


1 - _ 1 da 


1 f l dpb 


+ - ( —%-+3H)a 

Pb dt 


±^-(i + AL), ( 106 ) 

c 2 dt \ me 2 


A detailed study of the exact linearized equations 
(106)-(109) derived previously will be the subject of a 


future work. Here, we derive an approximate equation 
for the density contrast S(x, t) by neglecting the relativis¬ 
tic terms that involve a time derivative, i.e. terms like 
In this way, we simplify the equation for S(x,t) 
without altering the expression of the relativistic Jeans 
length. 

With this approximation, the linearized equations 
p6|-pl become 


dS 1- _ n 
— + -V- v = 0, 
dt a 


1 ( E 

~2 ( 1 4- 2 

c 2 V me 2 


h 2 


4 m 2 c 2 a 2 


AS- 1 


( 110 ) 


2c|A * 
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(111) 
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4m 2 c 2 y/pb 
3HH 2 dS 




4m 2 c 2 dt 


-a+H 


$ 


H 2 1 dpb < 9 $ 
m 2 c 4 pb dt dt 


, (107) 


h 2 


dv 

dt 
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4am 2 c 2 yfpb 


Hv=- h ^V(AS)--e 2 VS 

£(^ { )-4<v^) 


E 


rV<7 — 


amc* 


a \ c 


3 H 2 H , 
4am 2 c 2 
h 2 


„ , 1 dpb _<9$ 

V$ H-— V —, 

am 2 c 4 pb dt dt 


(108) 


dv 

dt 


n 2 - i 

Hv = — t V(AJ)-- 
4 m 2 a a a 


1 + 


2c 


V<F - -c 2 VJ, 
a 


( 112 ) 


A4> 

47t Gpbd 2 
(_ E 


= i + 


me 2 \2mc 2 


+ i U- 


6- 


E 


E 


me* 


+ 1] a 


me “ 
$ 

-n + 


E \ 2$ 
2 me 2 ) c 2 

3 H 2 ^ 

<f>. 


4-nGpbC 2 


(113) 


Taking the time derivative of Eq. (110) and the diver¬ 
gence of Eq. (112), we obtain 


c 
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2 H- 


dt 


= ~^A5 - 


h 2 


4m 2 a 4 


A 2 S 




(114) 
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Substituting Eq. ( |111| ) in Eq. ( |113| ), we get 

h 2 


. „ 2 — ( 1 + 2-| ) J — 

4n Gpba 2 ' 


+ 1 


$ 


E 


E 


me 2 V 2 me 2 


+ 1 


2c| 

C 2 J c 

2$ 

Ea 


E ( E 


2 me 2 \ 2 me 2 
$ 3 H 2 

-n + 


AirGpbC 2 


AS 

+ 1^6 

<1>. (115) 


Decomposing the perturbations in (spatial) Fourier 
modes, and using the identities (57) and (58) that re¬ 


main valid in the present context where now E(t) and 
c 2 (t) are time-dependent, we obtain a closed equation 
for Sk(t): 


AirGpb 

I 3 H 2 a 2 
k 2 c 2 


d 2 S k 

~W 


+ 2 H 


dS k 

dt 


h 2 k 4 


1 + 


2c 


1 + 


3c„ 


+ 


4m 2 a' 
h 2 k 2 
4m 2 c 2 a 2 


r 2 

4 +^k 2 

4 a z 


5k = 0 . 


(116) 


In a static universe (a = 1, H = 0) , taking Sk = 0 (i.e. 
w = 0), we can check that Eq. (116) reduces to Eq. (72) 
that determines the relativistic Jeans length. In the non- 
relativistic limit c -* +oo, the equation for the density 
contrast reduces to mm-- 


d 2 s k 

dt 2 


2 H 


dS k f K 2 k A c 2 
dt \4m 2 a 4 a 2 


k 2 — 47r Gpbj 5k = 0. 

(117) 


Because of the term appearing in the denominator of 
Eq. (116), the limits k —> 0 and c —> +oo do not com¬ 
mute. If we take the limit c —► +oo before the limit 
A: —>• 0, Eq. (116) reduces to 


g + 2/A)-4»C«J = 0 
dt z dt 


(118) 


which is equivalent to the CDM model where the evolu¬ 
tion of the perturbations is due only to the gravitational 
force. By contrast, if we take the limit k —> 0 before the 
limit c —> +oo, Eq. (116) reduces to 


=0 

dt 2 dt 


(119) 


in which the gravitational force has “disappeared”. This 
shows that the term appearing in the denominator of Eq. 
(116), which is purely relativistic, acts as an attenuation 


factor. Actually, this term can be written as 1 + (A/A h) 2 
where A = 27r /k is the wavelength of the perturbation 
and 


A h = 27t 


c 

y/3Ha 


( 120 ) 


is the cosmological horizon, or Hubble length. 10 This 
is the distance traveled by a photon with the velocity c 
during the Hubble time H -1 . The cosmological horizon 
sets the size of the observable universe. When A <C A h, 
which is the most relevant case, the term 1 + (A/A h) 2 
can be replaced by unity. However, when the wavelength 
A of the perturbation approaches the cosmological hori¬ 
zon, this term must be taken into account and it reduces 
the growth of the perturbations (it weakens the effective 
gravitational constant G e g(k) = G/{ 1 + k 2 H /k 2 )). As a 
result, it sets a natural upper cutoff, of the order of the 
Hubble length, for the perturbations that can grow. We 
shall come back to this reduction factor in Sec. IVII Cl 


VII. EVOLUTION OF THE DENSITY 
CONTRAST IN THE MATTER ERA 


In this section, we study the evolution of the den¬ 
sity contrast 5k in the matter era. We assume that the 
matter era starts at a t = 10 —4 , corresponding to the 
epoch of matter-radiation equality. This is just before 
the recombination epoch (matter-dominated era) start¬ 
ing at a rec = 1(V 3 when the electrons recombine with the 
atomic nuclei, after their last dispersion with the pho¬ 
tons, leaving an imprint of this interaction in the tem¬ 
perature fluctuations of the CMB. 


A. The equation for 5k (a) in the matter era 


In the matter era, we can m ake th e approximations 
~ PbC 2 and Pb e^, (see Sec. VC). In that limit, pb 


coincides with the rest-mass density and the Friedmann 
equations ([31]) , ([32]) and ([33]) reduce to 


^ + 3 H Pb = 0, 


( 121 ) 


H 2 = 


87 tG 


Pbi 


( 122 ) 


a 

a 


4t tG 


~Pb- 


(123) 


They lead to the EdS solution pi oc a 3 , a oc t 2 / 3 , H = 
2/(3 1) and pb = 1/(67 xGt 2 ). Using Eqs. (121)-(123), we 


can transform Eq. (116) for Skit) into an equation for 


10 Using the Friedmann equation ( |32| ), the Hubble length can be 
written as a sort of Jeans length A h = 2ire/ yj8nGeb/c^a in 
which the speed of sound is replaced by the speed of light. It can 
be interpreted as a maximum wavelength. 
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Sk{a) of the form 


B. Justification of the nonrelativistic limit 


d 2 S k 

da 2 

1 


3H 2 a 


k 2 c 


3 dS k 
2 a da 
, 2c 2 

c 2 


3 

2a? 


h 2 k 4 


T 


167rG/9f,m 2 a 4 47rGpf,a 2 


1 


3c 2 


h 2 k 2 


c 2 4m 2 c 2 a 2 


4 = 0, 


where we recall that c 2 is given by Eq. (591. Eq. (124) 


involves the following characteristic wavenumbers. The 
quantum Jeans wavenumber is defined by 


kn = 


{ 16 nGpbm 2 a' 


2„4\ !/ 4 


h 2 


(125) 


It can be written as kQ = Kqa 1 / 4 with Kq = 

(l6n G pba 3 m 2 /H 2 ) 1 ^ 4 . The self-interaction Jeans 

wavenumber is defined by 


It can be written as kj = Kja with kj = (Gm 3 / a s h 2 ^ 
The Compton wavenumber is defined by 


1/2 


kc = 


2 mca 
h 


(127) 


It can be written as kc = kc® with kc = 2 mc/H. The 
Hubble (or horizon) wavenumber is defined by 


kn = 


\/3Ha f 8nGpba 2 \ 




(128) 


It can be written as kn = kh /\fa with kh = 

1 /9 

(87t Gpba 3 /c 2 ) . We also introduce the dimensionless 
parameter 


4ira s h 2 pba 3 


(129) 


so that c 2 /c 2 = a/a 3 . In term of these parameters, Eq. 
(1241 can be rewritten as (we omit the subscript k on 
4]T 


(124) 


In this section, and in Secs. VII C|fVTI F we con¬ 
sider the case where the term associated with the Hubble 
length is very small: K 2 H /k 2 a <C 1. For a given wavenum¬ 
ber k, this corresponds to a scale factor a K 2 H /k 2 and, 
for a given scale factor a, this corresponds to a wavenum¬ 
ber k 3> Kn/a 4 / 2 . Therefore, this approximation is valid 
when the wavelength of the perturbation is much smaller 
than the cosmological horizon on the timescale consid¬ 
ered: A -C Atf(a). 

A direct (non-trivial) consequence of this assumption is 
that the term k 4 /Kqa associated with the quantum Jeans 
length is always much larger than the term k 2 / k 2 c a 2 as¬ 
sociated with the Compton length. Indeed, the condition 
k 4 /kqU k 2 / KqO , 2 is equivalent toa> Kq/k 2 K? c . Now, 


from Eqs. ( |125| ), (127) and (128) we can easily establish 
the identities 

( k c k H ) 2 = 2 kq and ( k c k h ) 2 = 2 Kq. (131) 

As a result, Kq/k 2 Kg. = k 2 h / 2k 2 < K 2 H /k 2 . There¬ 
fore, the condition a K 2 H /k 2 automatically implies 
a Kq/k 2 K,Q which itself implies k 4 1 Kqa k 2 / KqO, 2 . 
We stress that this inequality is independent on the mass 
m of the bosons. 

On the other hand, the condition a/a 3 <C 1 will be 
satisfied during the whole matter era if it is satisfied ini¬ 
tially, i.e. if a <C a 3 = 10 -12 . This condition corresponds 
to c s <C c. In Appendix [G] we show that this condition is 
always satisfied in the matter era. 

In conclusion, if we assume that the wavelength of the 
perturbation is always much smaller than the cosmolog¬ 
ical horizon during the matter era ( A <C A //(a)), and if 
we use the fact that a/a 3 <C 1, Eq. (130) reduces to 


d 2 d 3 dS 
da 2 + 2a da + 2a 2 


AW 747 + 70 —1H = 0, (132) 


K 4 Q a 


K 2 7 a 2 


which is the nonrelativistic equation for the density con¬ 
trast of a SF/BEC |97j, 100] . Therefore, the nonrelativis¬ 
tic limit is valid in the matter era when A <C A h(o). 

The CDM model is recovered by taking Kq — > +oo and 
kj —^ Too in Eq. (132) yielding 


d 2 S, 


CDM 


da 2 


3 dScDM 3 

2cl ~da 2a?^ CDM = °’ 


(133) 


<PS_ 
da 2 


3_dS _3_T_fc^ fc 2 

2 a da 2a 2 [Kqa K 2 a 2 


IT 


k 2 a 




+ 


k 2 A 

4° 2 /. 


(5 = 0. 
(130) 


It is interesting to see how the different characteristic 
scales Kq , kj, kjj and kc naturally appear in this equa¬ 
tion. 


The growing solution is Sc dm oc a. It is usually con¬ 
sidered that Si ~ 10 -5 at the initial time a* = 10~ 4 
of matter-radiation equality [3]. Therefore, the growing 
evolution of the density contrast in the CDM model can 
be written as 

c 

Scdm{o) = —a. (134) 

a I 

We will take this CDM result as a reference, and compare 
it with the prediction of the SF model. We note that 
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Eq. ( 132 ) for the density contrast of the SF reduces to 
Eq. (1331 when k —» 0 and when a —» +oo because 
the quantum term k 4 / Kqa and the self-interaction term 


k 2 /k 2 jCL 2 become negligible in these limits. Therefore, the 
SF is expected to behave similarly to CDM at large scales 
and at late times. 


C. Non-interacting scalar field 


For a noninteracting SF (kj —> +oo), the equation 
determining the evolution of the density contrast reduces 
to 



- 1 \S = 0 . 


(135) 


It turns out that Eq. (135) can be solved analytically 


m- Its growing solution is given by 


6(a) = —tjtJ 6 
v ’ a 1 / 4 


V6- 


k q a 


1/2 


(136) 


We shall determine the amplitude A so that the asymp¬ 
totic behavior of the density contrast of the SF at late 
times exactly matches the solution (134) of the CDM 


model. Using the asymptotic expansion of the Bessel 
function for small arguments (see, e.g., Appendix A of 
)100p . we obtain 


A = 


6 5 / 4 [tt Si f k 
2 a , ; \ kq 


(137) 


Eq. (136) with A given by Eq. (137) determines the 


evolution of the density contrast of a non-interacting SF 
during the matter era in the nonrelativistic limit A -C 
A if (a). 

We can identify two regimes. We first consider the 
case k 4 / k 4 qU ^ 1- For a given wavenumber k , this cor¬ 
responds to a scale factor a <C k 4 /K,Q. Alternatively, for 
a given scale factor a , this corresponds to a wavelength 
A -C Aq(ci). Since the wavelength of the perturbation 
is smaller that the quantum Jeans length, the density 
contrast 5(a) displays oscillations. Using the asymptotic 
expansion of the Bessel function for small arguments, we 
find that 


a fc 4 /«Q, the perturbation behaves like in the CDM 
model, see Eq. (134). 


In summary, when a -C k 4 / Kq, i.e. A <C Ag(a), the 
perturbation oscillates with a constant amplitude and 
when a k 4 / Kq , i.e. A A q(o), the perturbation 
grows linearly with the scale factor like in the CDM 
model. 


To fix the ideas, we consider that the matter era starts 
at Oj = 10~ 4 and we study the evolution of the density 
contrast up to ao = 1, corresponding to the present-day 
universe (of course, our linear analysis is not valid up to 
that point, but our aim is just to give a simple illustra¬ 
tion of our results). Our nonrelativistic treatment is valid 
on this interval provided that m ZS> ti 2 H /k 2 . This corre- 
sponds to wavenumbers k kh/ a/ . This constraint is 
not restrictive since Kh/ aj 2 = 6.49 x 10 - 25 m -1 is very 
small (see Table I in Appendix [f|. For k/ kq > a ^ 4 = 1, 
the evolution of the density contrast is purely oscilla¬ 
tory (see Fig. [i). For k/nQ < a. 4,/4 = 0.1, there is no 
oscillation and the density contrast grows from the be¬ 
ginning (see Fig. [3]). 11 In the intermediate situation 
a 4 / 4 = 0.1 < k/ kq < a 4/4 = 1, the density contrast first 
presents oscillations until a ~ k 4 / Kq, then grows (see 
Fig. 0. We note that the density contrast in the SF 
model grows faster than in the CDM model. 



FIG. 2: Evolution of the density contrast 5(a) of a non¬ 
interacting SF for k/KQ = 1 (semi-log plot). It is compared 
with the density contrast 5c dm (a) of the CDM model. 


6(a) ~ -2— f — 'j sin (y /&—2 1721 • ( 138 ) 

a i \ K Q / \ k q a ' J 

This solution is valid for a k 4 / Kq. We see that the am¬ 
plitude of the oscillations is constant, i.e. independent on 
the scale factor. We now consider the case k 4 /KqO 1. 
For a given wavenumber k , this corresponds to a scale 
factor a k 4 / Kq. Alternatively, for a given scale factor 
a , this corresponds to a wavelength A Aq(ci). Since the 
wavelength of the perturbation is larger that the quan¬ 
tum Jeans length, the density contrast 5(a) increases. For 


The SF has several virtues. First, it does not create an 
over-abundance of small-scale structures, contrary to the 
CDM model, because the perturbations below the quan- 


11 We can give a more general criterion for the absence of oscil¬ 
lations, independent on the interval of time considered. Our 
nonrelativistic treatment requires that a /k 2 and the den¬ 

sity contrast displays oscillations provided that a <g fc 4 /«g. 
Therefore, the density contrast will not display oscillations if 

































































17 



FIG. 3: Evolution of the density contrast 5(a) of a non¬ 
interacting SF for k / kq = 0.1 (semi-log plot). 



a 


This equation can be solved analytically m • The grow¬ 
ing solution is given by 


^( a ) ~ a i/4^-| 


(M) 


(140) 


The amplitude B is determined by requiring that the 
asymptotic behavior of Eq. (140) for a —>• +oo exactly 


matches the solution (1341 of the CDM model. This gives 


b = r — 


3 5/S 

4 ; 2 15 /s 


-) 

KjJ 


5/4 


(141) 


Eq. (140) with B given by Eq. (141) determines the 


evolution of the density contrast of a self-interacting SF 
during the matter era in the nonrelativistic limit A <C 
A h{o) and in the TF limit. 

We can identify two regimes. We first consider the 
case k/nja 1. For a given wavenumber k, this cor¬ 
responds to a scale factor a <C k/Kj. Alternatively, for 
a given scale factor a, this corresponds to a wavelength 
A -C Aj(a). Since the wavelength of the perturbation is 
smaller that the self-interaction Jeans length, the density 
contrast 6(a) displays oscillations. Using the asymptotic 
expansion of the Bessel function for small arguments, we 
find that 


W~r(-i) 


3 3 / 8 1 6 t 1/4 / k A 3/4 

2 9 / 8 ypK oq a \ Kj J 


x cos 


/ f3k_l 

\ \ 2 kj a 



(142) 


FIG. 4: Evolution of the density contrast 5(a) of a non¬ 
interacting SF for k/ttQ = 0.5 (semi-log plot). 


turn Jeans length A q do not grow (they oscillate). The 
quantum Jeans length provides therefore a sharp small- 
scale cut-off in the matter power spectrum. On the other 
hand, the perturbations above the quantum Jeans length 
Aq grow sensibly more rapidly than in the CDM model. 
This can accelerate the formation of large-scale struc¬ 
tures in the nonlinear regime in agreement with some 
observations m- 


D. Self-interacting scalar field in the 
Thomas-Fermi limit 


We now consider a self-interacting SF. We make the 
Thomas-Fermi (TF) approximation (kq —> +oo) which 
amounts to neglecting the quantum term in front of the 
self-interaction term. In that case, the equation deter¬ 
mining the evolution of the density contrast reduces to 


<P8_ 3 / k 2 \ 

da 2 + 2a da~^~ 2a 2 A« 2 a 2 ) 


(139) 


This solution is valid for a -C k/nj. We see that the 
amplitude of the oscillations grows as a 1 / 4 . This is a dif¬ 
ference with respect to the non-interacting SF for which 
the amplitude of the oscillations is constant (see Sec. 


VIIC). We now consider the case k/Kja <C 1. For a 


given wavenumber k , this corresponds to a scale fac¬ 
tor a k/nj. Alternatively, for a given scale factor a, 
this corresponds to a wavelength A /?> A j(a). Since the 
wavelength of the perturbation is larger that the self- 
interaction Jeans length, the density contrast <5(a) in¬ 
creases. For a /§> k/ kj, the perturbation behaves like in 
the CDM model, see Eq. (134). 

In summary, when a <C k/ kj, i.e. A <C A j(a), the per¬ 
turbation oscillates with a growing amplitude scaling as 
a 1,/4 and when a k/nj , i.e. A A j(a), the perturba¬ 
tion grows linearly with the scale factor like in the CDM 
model. 

As before, we consider that the matter era starts at 
ai = 10” 4 and we study the evolution of the density 
contrast up to ag = 1. Our nonrelativistic treatment is 
valid for k 3> Kn/a^ 2 = 6.49 x 10 _25 m _1 . For k/nj > 
cio = 1, the evolution of the density contrast is purely 
oscillatory (see Fig. [5]). For k/nj < ai = 10 -4 , there 
is no oscillation and the density contrast grows from the 
beginning (the evolution is similar to the one represented 
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in Fig. [3]). 12 In the intermediate situation a, = 10” 4 < 
h/nj < 1, the density contrast first presents growing 
oscillations until a ~ k/nj, then grows linearly with a 
(see Fig. [6]). We note that the density contrast in the SF 
model grows faster than in the CDM model. 


cd 

CO 



FIG. 5: Evolution of the density contrast 8(a) of a self¬ 
interacting SF in the TF limit for k/nj = 1 (semi-log plot). 



FIG. 6: Evolution of the density contrast 8(a) of a self¬ 
interacting SF in the TF limit for k/nj = 0.1 (semi-log plot). 


length A j. This is an advantage of the SF model with re¬ 
spect to the CDM model. Note that the self-interaction 
Jeans length, which provides a sharp small-scale cut-off 
in the matter power spectrum, corresponds precisely to 
the size of the dark matter halos (see Eq. (G2|) that 
form in the nonlinear regime [5Z1I67]. 

Remark: we note that the Compton term in Eq. (130) 
scales as —k 2 /K^a 2 similarly to the s elf-interaction Jeans 
term k 2 /n 2 a 2 . Therefore, Eq. (139) with the Compton 
term retained can still be solved analytically by simply 
replacin g 1 / k 2 by 1/k 2 — 1 /k (G How ever, we have shown 
in Sec. |VII B (see also Sec. VUG) that the Compton 
term is always negligible in the matter era. 


E. Self-interacting scalar field beyond the TF limit 


We still consider a self-interacting SF but we now take 
the quantum term into account. The evolution of the 
density contrast is determined by Eq. (132). This equa¬ 


tion has to be solved numerically. We would like to com¬ 
pare the solution of Eq. (1321 with the solution ( 134| ) of 
the CDM model by imposing that they coincide at late 
times. However, there is a difficulty. Since the solution 
of Eq. (132) may be strongly oscillating at early time, 


it is not convenient to specify the initial condition at 
and solve the equation forward in time. Such a procedure 
would be very sensitive to the initial condition and would 
not in general reproduce the CDM result at late times. 13 
Instead, we can specify the initial conditions at ao = 1 
(today) and solve the differential equation backwards in 
time until at = ICC 4 . This procedure has already been 
used in the past [55j. Therefore, we solve Eq. (132) 


numerically with the initial condition (So, So) obtained 
from the analytical solution determined by Eqs. (136) 
and ( |137 ) at ao = 1. For illustration, we have taken 
k/nQ = 0.5 and kq/kj = 0.1 in Fig. [?] 

The evolution of the perturbation can be understood 
qualitatively by using simple scaling arguments. The 
self-i nteraction term dominates the quantum term in Eq. 
(132) when k 2 /n 2 a 2 k 4 /KqCL- For a given wavenum¬ 
ber k, this corresponds to a scale factor a <C Hq/^jk 2 . 
Alternatively, for a given scale factor a, this corresponds 


The remarks concerning the virtues of the non¬ 
interacting SF made at the end of Sec. 


VIIC 


can 

be transposed to the self-interacting SF with the self¬ 
interaction Jeans length replacing the quantum Jeans 
length. The SF suppresses small-scale structures below 
the self-interaction Jeans length A j and accelerates the 
formation of structures above the self-interaction Jeans 


12 The more general criterion for the absence of oscillations, similar 

1/3 2/3 

to the one obtained in footnote 11, is k <C kJ . 


13 This is a particularity of the SF. Because of the oscillations of 
<5(a) at early time, the evolution of the perturbations with small 
wavelengths A is very dependent on the initial values of ((5^,(5*) 
prescribed at a^. We need very particular values of (<5^, <5*) in 
order to recover the results of the CDM model at later times. 
This implies that, in general, the evolution of the perturbations 
with small wavelengths is different in the SF and CDM models, 
even at late times. This is not a drawback of the SF model 
because we know that the CDM model presents discrepencies 
with observations at small scales and that the SF model has 
been introduced precisely to solve them (see the discussion at 
the end of Secs. | VII C| and | VII D| . Of course, for sufficiently 
large wavelengths A, the fine tuning problem disappears, and the 
SF and CDM models give the same results at late time. 
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to a wavelength A Xq/Xj. For a <C Kq/kj/c 2 , the 
system is in the TF regime and for a Kq/ n 2 jk 2 it is in 
the non-interacting regime. In the TF regime, the per¬ 
turbation undergoes growing oscillations for a <C k/nj 
and grows linearly for a k/ kj. In the non-interacting 
regime, the perturbation undergoes oscillations with a 
constant amplitude for a -C k 4 / Kq and grows linearly 
for a k 4 /k/q. 



ing oscillations (provided that fc/rej > dj) followed by 
a linear growth. The growing oscillations take place at 
the begining of the TF regime. The linear growth starts 
at the end of the TF regime and continues in the non¬ 
interacting regime (there is no oscillation in the non¬ 
interacting regime because k 4 / Kq < Kq/K 2 k 2 ). It is not 
possible to clearly see the change of regime during the 
linear growth so the evolution of the perturbation looks 
similar to that reported in Fig. [6] 

(ii) If fc 4 / Kq > K^g/K^jk 2 , we necessarily have k/nj > 
Kq/Kjk 2 . In that case, the perturbations undergo grow¬ 
ing oscillations, followed by constant-amplitude oscil¬ 
lations, and finally by a linear growth (provided that 
k 4 /Kq < 1). The growing oscillations take place in the 
TF regime. The constant-amplitude oscillations and the 
linear growth take place in the noninteracting regime 
(there is no linear growth in the TF regime because 
k/Kj > Kq/k 2 /c 2 ). In that case, the change of regime is 
clearly apparent since we can see the difference between 
growing oscillations and constant-amplitude oscillations. 
This is illustrated in Fig. [8] 

F. The case of a negative scattering length 


FIG. 7: Evolution of the density contrast 5(a) of a self¬ 
interacting SF for k/KQ = 0.5 and kq/kj = 0.1 (semi-log 
plot). It is compared with the non-interacting SF solution 
(dashed line). 



FIG. 8: Evolution of the density contrast 5(a) of a self¬ 
interacting SF for k/ kq = 0.840 and kq/kj = 0.022 (semi-log 
plot). 

If Kq / Kjk ' 2 <C a t . the system is always in the non¬ 
interacting regime. If Kg/«; 2 fc 2 clq = 1, the system is 
always in the TF regime. In these cases, we are led back 
to the situations considered in Secs. IVII Cl and IVII Dl 
We now assume a* < Kq/Kjk 2 < 1. Two situations can 
occur. 

(i) If k/nj < KQ/n 2 jk 2 , we necessarily have k 4 /Kq < 
Kq/ n 2 jk : 2 . Therefore, the perturbation undergoes grow- 


The previous results implicitly assume that the scat¬ 
tering length is positive (a s > 0), corresponding to a 
repulsive self-interaction. In this section, we consider the 
case of a negative scattering length (a s < 0), correspond¬ 
ing to an attractive self-interaction. The equation for the 
density contrast becomes 


cPS 3 d5 

da 2 ^ 2a da ^ 2 a 2 \ K“ka kPjO 2 


k 4 

~A 

V Q 


AW -TT-WW-1 U = 0, (143) 


where now kj = (Gm 3 /|a s |fi. 2 ) 4 / 2 . We first neglect the 
quantum term in Eq. ( |143[ ) (TF approximation). In that 
case Eq. (1431 reduces to 


d H 

da 2 


~ “-b 7TV7 (-2 2 — 1 ) — 0. 

Kj a 


3 dS 3 
2a da 2o 2 


The growing solution of Eq. (144) is given by 
1 


S(a) oc -^I_ 


3 k 1 
2 Kj a 


(144) 


(145) 


This solution diverges exponentially rapidly towards —oo 
when a 0. This divergence can be avoided by tak¬ 
ing the decayi ng so lution into account TOOL Therefore, 
we solve Eq. (144) with the initial condition Si = 10 -5 
and Sj = Si/ai = 0.1 at a* = 10 —4 , obtained from Eq. 
(1341. The result is represented in Fig. [9] For small 
k (large scales), the density contrast evolves as in the 
CDM model. For large k (small scales), the density con¬ 
trast increases exponentially rapidly with the scale factor 
because the attraction due to the self-interaction adds 
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to the gravitational attraction [lOOj . An attractive self¬ 
interaction therefore favors the growth of structures, but 
it does it in a rather dramatic manner. Indeed, it leads 
to an over-abundance of substructures at a very early 
epoch. This is in contradiction with the observations 
that do not show such an abundance of substructures. 
The CDM model, in which only gravity is in action, al¬ 
ready presents too many substructures. Therefore, we 
need a mechanism that erases small-scale structures, not 
a mechanism that enhances them. Consequently, the SF 
should have in principle a positive scattering length, not 
a negative scattering length. 14 



a 


FIG. 9: Evolution of the density contrast 5(a) of a self¬ 
interacting SF with a s < 0 in the TF limit for (bottom to 
top) k/nj = 0, k/nj = 5 x 10~ 4 and k/nj = 10 -3 (log-log 
plot). 


These conclusions are not true anymore if |a s | is ex¬ 
tremely small (typically |a s | -C 10 -60 fm) because, in 
that case, we need to take into account the quantum term 
in Eq. (1431 that stabilizes the structures at small scales. 
The quantum term becomes efficient when a k'q/K 2 jk 2 . 
Therefore, it rapidly stabilizes the modes with large k 
(small scales) and stops their growth. Indeed, after an 
exponential growth, the perturbation starts to oscillate 
(see Fig. [To] ) . Modes with intermediate k (intermedi¬ 
ate scales) grow exponentially rapidly because the range 
where the attractive self-interaction dominates is longer 
(i.e. Kq/Kjk 2 is larger). Finally, modes with small k 
(large scales) grow slower (linearly with a) and behave 
like the CDM model. Therefore, bosons with an ex¬ 
tremely small negative scattering length may be interest¬ 
ing DM candidates (as suggested in jlOOj ) because they 
can accelerate the growth of perturbations, hence the 
formation of structures, at intermediate scales without 


14 In addition, it is shown in 1571 1581 that it is not possible to 
form realistic dark matter halos when a s < 0 because they are 
unstable above a maximum mass M max = 1.01 '2hj y 7 Grn \a s that 
is usually very small. 


forming too many substructures at small scales. These 
considerations may be relevant to QCD axions because 
they usually have an extremely small negative scattering 
length [521 . 



FIG. 10: Evolution of the density contrast 5(a) of a self¬ 
interacting SF with a a < 0 for e = kq/kj = 0.0085 and 
k/kQ = 0.7 (semi-log plot). The dashed curve (exponen¬ 
tial growth) represents the solution of Eq. ( |143| ) without the 
quantum term and the solid curve shows the stabilization due 
to the quantum term. 


G. Relativistic attenuation when A approaches the 
cosmological horizon 

We now consider the case where the wavelength A of 
the perturbation approaches the cosmological horizon A h 
so that relativistic effects must be taken into account. 

If we consider wavelengths A ~ Aj y, or larger, then 
the Compton term k 2 /n 2 c a 2 in Eq. (1301 is negligible 
in front of unity (k 2 /n^a 2 <C 1, i.e. A A c) because 
A h A c (see Appendix 0- Since we have shown in 
Sec. | VII B| that the Compton term is negligible in front 
of the quantum term when A -C A h, we conclude that 
the Compton term is always n egligibl e in the matter era. 
We have also shown in Sec. VIIB that a/a 3 -C 1 in 


the matter era. Therefore, the evolution of the density 
contrast in the matter era, valid for any perturbation, is 
determined by the equation 


d 2 5 3 dS 3 
da 2 + 2a da 2a 2 


fc 4 

K Q a 


k 2 


1 n 1 


k , a 


1 

A 

1 + vF 

k 2 a 


5 = 0. 


(146) 


In the matter era, we have A c (A q, A j) <C A h (see 
Appendix [F]). For A -C (Aq,Aj), we can use the nonrel- 
ativistic results of Secs. |VII C|fVtl E| since A -C A h- In 
that case, we have seen that the perturbation oscillates 
because the wavelength is smaller than the Jeans length. 
We now assume A (Aq,Aj). In that case, Eq. (146) 
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FIG. 11: Evolution of the density contrast <5(a) of a self¬ 
interacting SF for k/nn = 10 and k/KH = 0.1 (log-log 
plot). We have also plotted the CDM soluti on ( recovered 
for k kh / \[a) and the asymptotic solution (1491 valid for 
k ->• 0. 


can be approximated by 


cP8_ 

da 2 


3 dd 
2 a da 


2 a 2 


-6 = 0 . 


1 + 


(147) 


k 2 G 


We solve Eq. (1471 with the initial condition 8i = 10 5 
and Si = 6i/ai = 0.1 at cq = 10 -4 , obtained from Eq. 


(1341. The last term in Eq. (147) corresponds to the 


gravitational attraction. For A <C A h, we can use the 
results of Secs. |VIIC|VilE| and we find that the per¬ 
turbation grows linearly with a, exactly like in the CDM 
model . In deed, the term n 2 H /k 2 a can be neglected and 
Eq. ( |147[ ) reduces t o Eq . (133) which has the grow¬ 
ing solution of Eq. (|134[) 


We now consider the situa¬ 
tion where A approaches the Hubble length A h- la that 
case, the weight of the gravitational term is reduced be¬ 
cause of relativistic effects and, consequently, the pertur¬ 
bations grow slower. This is illustrated in Fig. mi When 
A —I +oo , the gravitational term becomes negligible and 
Eq. (11471) reduces to 


d 2 S 

da 2 


~= 0 . 
2 a da 


It has the solution 


6(a) =6i ( 3-2J-A 


(148) 


(149) 


which tends to a constant 3<5i when a —> +oo. This 
shows that when the wavelength approaches or overcomes 
the cosmological horizon, the perturbations do not grow 
anymore. There is a attenuation of relativistic origin. 

Coming back to the general equation (146) and sum¬ 


marizing our results, we have shown that the relativistic 
term l + k^/k 2 sets a natural upper cutoff, of the order 
of the Hubble length A h, above which the perturbations 


are attenuated. On the other hand, the terms k 4 /kq and 
k 2 /kq set natural lower cutoffs, of the order of the Jeans 
length A q or A j, below which the perturbations oscillate. 
In between, the perturbations grow linearly with a as in 
the CDM model. Therefore, the effect of the relativis¬ 
tic SF is to introduce sharp cut-offs in the matter power 
spectrum at small scales A q or A j because of quantum 
mechanics and at large scales A# because of general rel¬ 
ativity. 


VIII. ABOUT THE IMPORTANCE OF 
RELATIVISTIC EFFECTS 

Before concluding this paper, we would like to discuss 
the importance, or non-importance, of relativistic effects 
in the SFDM model. 

In the matter era, we have shown that relativistic ef¬ 
fects can be neglected as long as the wavelength of the 
perturbations is much smaller than the cosmological hori¬ 
zon. This is the case for most perturbations of physical 
interest. Since there is a wide range of scales satisfying 
(Aq,Aj) -C A <C A h, the perturbations can grow during 
the matter era and lead to the formation of structures 
in the nonlinear regime. However, in the stiff matter era 
and in the radiation era, relativistic effects are impor¬ 
tant because the wavelength of the perturbations may 
more easily reach the horizon (which is smaller). Since 
general relativity tends to prevent the growth of pertur¬ 
bations when their wavelength approaches the horizon, 
this may explain why the perturbations do not grow in 
the radiation era (the Jeans length may be larger than 
the horizon). This interesting effect will be studied in 
detail in future works. 

On the other hand, even if relativistic corrections are 
weak in the matter era, we may wonder whether their ef¬ 
fect could be detected in an era of “precision cosmology” 
11201 . In particular, it would be interesting to see if one 
can observe differences between the KGE equations con¬ 
sidered in this paper and the heuristic KGP equations 
studied in the past, in which gravity is introduced by 
hand in the KG equations (see Appendix [d|) . 

Relativistic effects are important for other astrophysi- 
cal systems described by SFs or BECs. We can mention, 
for example, the case of boson stars IHH32 and the case of 
microscopic quantum black holes made of BECs of gravi¬ 
tons stuck at a quantum critical point (szniiinng. it 
has also been proposed |32] that, because of their super¬ 
fluid core, neutron stars could be BEG stars. Indeed, 
the neutrons (fermions) could form Cooper pairs and be¬ 
have as bosons of mass 2 m n . This idea may solve cer¬ 
tain issues regarding the maximum mass of neutron stars. 
To be complete, we should also mention analog models 
of gravity in which condensed matter systems described 
by the GP equation or by the KG equation are used to 
simulate results of classical and quantum field theory in 
curved space-time (12311125] . One of the first models for 
experimentally simulating black hole (Hawking) evapora- 




















22 


tion was suggested by Unruh in 1981 jl25j . Fifteen years 
later, a sonic analog of black holes specifically using BECs 
was developed 11261 . Since then, various configurations 
of relativistic BECs have been proposed to simulate dif¬ 
ferent scenarios of gravity |127H129j . 


IX. CONCLUSION 

In this paper, we have developed a formalism based on 
a relativistic SF described by the KGE equations in the 
weak field limit. We have transformed these equations 
into equivalent hydrodynamic equations. These equa¬ 
tions are arguably more tractable than the KGE equa¬ 
tions themselves. In the nonrelativistic limit, they re¬ 
duce to the hydrodynamic equations directly obtained 
from the GPP equations m■ Therefore, our formal¬ 
ism clarifies the connection between the relativistic and 
nonrelativistic treatments. We note that, in the rela¬ 
tivistic regime, the hydrodynamic variables ip, p, S, v, 
p,... that we have introduced do not have a direct phys¬ 
ical interpretation. It is only in the nonrelativistic limit 
c —> Too that they coincide with the wavefunction, rest- 
mass density, action, velocity, and pressure. However, 
these variables are perfectly well-defined mathematically 
from the SF ip in any regime, and they are totally legiti¬ 
mate. Furthermore, in terms of these variables, the rela¬ 
tivistic hydrodynamic equations take a relatively simple 
form that provides a natural generalization of the nonrel¬ 
ativistic hydrodynamic equations. The main advantage 
in using the hydrodynamic representation of the SF is the 
facility of interpreting the results since the hydrodynamic 
equations are expressed in terms of familiar physical vari¬ 
ables (density, velocity, pressure...). In addition, the use 
of these variables avoids the need of making averages over 
the oscillations of the SF as discussed in Sec. |V] This is 
valuable because it is not easy in practice to deal with 
the oscillations of the SF when directly solving the KGE 
equations numerically. 

The complete study of the relativistic hydrodynamic 
equations is of considerable interest but it is, of course, 
of great complexity. In this paper, we have started their 
study in simple cases and we have obtained explicit re¬ 
sults. Our paper shows therefore that these equations 
can be useful in practice. 

First, we have considered the relativistic Jeans insta¬ 
bility problem in a static universe. We have derived the 
relativistic dispersion relation and the relativistic Jeans 
length. Physical applications of these results will be de¬ 
veloped in our companion papers. 

Secondly, we have checked that the hydrodynamic 
equations of the SFDM model reproduce the evolution of 
the homogeneous background obtained previously by Li 
et al. lfl 1 directly from the KGE equations: a stiff mat¬ 
ter era, followed by a radiation era (for a self-interacting 
SF), and a matter era. 

Finally, we have started to analyze the evolution of 
the perturbations in an expanding universe dominated 


by a SF. We have considered the linear regime and we 
have focused on the matter era during which the growth 
of structures is expected to take place. In this regime, 
we can use the nonrelativistic limit of our formalism to 
treat the perturbations that are well inside the cosmolog¬ 
ical horizon (A <C A h)- We have evidenced analytically 
and numerically different regimes. Perturbations with 
A < Aj oscillate with an amplitude growing as a 1 / 4 , 
perturbations with A <C A q oscillate with a constant am¬ 
plitude, perturbations with (Aq,Aj) < A « A h grow 
linearly with a like in the CDM model, and perturba¬ 
tions with A ~ A h are attenuated by general relativity 
and do not grow anymore. Therefore, the relativistic SF 
introduces sharp cut-offs in the matter power spectrum 
at small scales (Aq or A j) because of quantum mechanics 
and at large scales (A h) because of general relativity. 

Generically, we can use the TF approximation at early 
times, then the non-interacting approximation, and fi¬ 
nally the classical (non quantum) approximation. As 
a result, the perturbations first display growing os¬ 
cillations, followed by constant amplitude oscillations, 
and finally grow linearly with a similarly to the CDM 
model (see Fig. [8]). The growth of perturbations in 
the SF model is substantially faster than in the CDM 
model in agreement with certain cosmological observa¬ 
tions where large-scale structures are observed at high 
redshifts m- Depending on the wavelength of the per¬ 
turbations, the duration of these different regimes in the 
matter-dominated era can change. We have provided all 
the scaling laws necessary to have a complete picture of 
the evolution of the perturbations in any case. In future 
works, we will extend our study of perturbations in the 
radiative era and in the stiff matter era where relativistic 
effects are expected to play an important role. 

In this paper, we have studied the development of per¬ 
turbations of the SFDM model in the linear regime. In 
future works, it will be important to consider the non¬ 
linear regime where structure formation actually occurs. 
In general, this problem must be addressed numerically. 
The hydrodynamic equations derived in this paper may 
be very helpful because they may be easier to solve than 
the KGE equations. 15 Therefore, numerical simulations 
using fluid dynamics should be developed in the future. 
As a first step, relativistic effects could be neglected and 
the nonrelativistic equations of Appendix [C] could be 
used. 16 These equations are similar to the hydrodynamic 
equations of CDM except that they include a quantum 


15 Since our equations are based on the Newtonian metric assum¬ 
ing that <T>/c 2 <C 1, they cannot be extrapolated to the nonlinear 
relativistic regime where terms like <1> 2 /c 4 must be taken into ac¬ 
count. However, our equations are exact in the linear relativistic 
regime which is relevant to describe the radiation era, and in the 
nonlinear nonrelativistic regime (c —> Too) which is relevant to 
study the formation of structures in the matter era. 

16 It may be mentioned that very nice numerical simulations 
have been made recently for nonrelativistic non-interacting 
SFDM/BECDM by directly solving the SP equations 16011611 ■ 
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potential (Heisenberg) and a pressure term (scattering) 
that avoid singularities at small scales and may solve the 
missing satellite problem [lOOj . In this respect, it may 
be recalled that the SP equations were introduced early 
by Widrow and Kaiser USD] as a procedure of small-scale 
regularization (a sort of mathematical trick) to prevent 
singularities in collisionless simulations of classical parti¬ 
cles. In their approach, H is not the Planck constant, but 
rather an ajustable parameter that controls the spatial 
resolution. Their procedure finds a physical justification 
if DM is made of self-gravitating BECs |131] . 

In this paper, we have studied the evolution of a uni¬ 
verse containing exclusively a SF representing dark mat¬ 
ter. If we want to construct a more general cosmologi¬ 
cal model, we must take into account the contribution of 
other components such as normal 17 radiation (due to rel¬ 
ativistic particles: photons, neutrinos...), baryonic mat¬ 
ter, dark energy (e.g. the cosmological constant), and 
possibly non-gravitational couplings (e.g. electromag¬ 
netism). Normal radiation, baryons, and the cosmologi¬ 
cal constant can be introduced straightforwardly in the 
energy-momentum tensor or in the Friedmann equations. 
Since our perturbation analysis is carried out during the 
matter-dominated era (dark matter), corresponding to 
10~ 4 < a < 1, the contribution to the density contrast 
of the components different from the SF is subdomi¬ 
nant, so these components are not expected to change 
the evolution of the perturbations significantly. In the 
radiation-dominated era, corresponding to a < 10~ 4 , the 
growth of structures is inhibited because the velocity dis¬ 
persion of the particles dominates their gravitational at¬ 
traction (free streaming) and the wavelength of the per¬ 
turbations that could trigger gravitational collapse lies 
outside the horizon making these perturbations ineffec¬ 
tive. Non-trivial gravitational couplings must be treated 
specifically. For example, electromagnetism can be taken 
into account by adding the contribution of the Maxwell 
tensor in the Lagrangian. The hydrodynamic equations 
can then be generalized. In the case of a complex SF 
(charged), electromagnetic forces can have a significant 
influence on the formation of structures due to vortical 
motion. 

We would like to conclude this paper by recalling the 
interest of the SFDM model with respect to the CDM 
model for what concerns astrophysical and cosmological 
observations. One of the main observations that can be 
related to our study is the recently observed large-scale 
structures at high redshifts [119j . Contrary to the CDM 
model which does not seem to give sufficient time for 
such structures to form, SFDM perturbations seem to 
grow somewhat faster, thereby facilitating the formation 
of structures at earlier times, or more rapidly. Another 


17 We call it “normal” in order to distinguish it from the radia¬ 
tion produced by a relativistic self-interac ting S F, which has a 
different physical origin (see 11041 and Sec. VC I. 


observation that can be related to our work is the “miss¬ 
ing satellite problem” m- Numerical simulations of 
CDM lead to an over-abundance of small structures be¬ 
cause the Jeans length is equal to zero (or is very small). 
Such a large number of satellite galaxies is not observed 
in the Universe. By contrast, the SFDM model provides 
a natural cutoff for the formation of small structures, 
due to quantum mechanics, giving a possible solution to 
the missing satellite problem because the Jeans length is 
finite. Another problem of the CDM model is the “cusp- 
core problem”. The CDM model predicts that DM halos 
should be cuspy (the central density should diverge as 
r _1 ) il33j while observations reveal that they have a flat 
core 1134!. Once again, the quantum properties of the 
SFDM model may solve this problem. Indeed, the quan¬ 
tum potential arising from the Heisenberg uncertainty 
principle (for a noninteracting SF) or the pressure due to 
the scattering (for a self-interacting SF) prevent gravita¬ 
tional collapse at small scales and lead to central density 
cores instead of cusps. Other astrophysical and cosmo¬ 
logical applications of our SFDM formalism will be given 
in future papers. 


Appendix A: The value of A 


The GP equation is obtained from the KG equation by 
means of the transformation 


ip = Ae 


2< /V- 


(Al) 


The constant A can be computed as follows. Substituting 
Eq. (Al) in Eq. (21), we find that the energy density of 
the SF is given by 
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Taking the nonrelativistic limit c —> +00 of this equation, 
we obtain 


e m 2 A 2 2 m 2 A 2 


(A3) 


where p = \ip\ 2 is the rest-mass density. Since e ~ pc 2 in 
the nonrelativistic limit c -A +00, we find 



(A4) 


Appendix B: Spatially inhomogeneous clusters 

We consider a spatially inhomogeneous SF/BEC clus¬ 
ter at equilibrium, representing for example a dark mat¬ 
ter star, a boson star, or a dark matter halo. This cluster 
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results from the nonlinear development of the Jeans in¬ 
stability. It corresponds to a steady state of the hydrody¬ 
namic equations < 41)-(44) with a = 1, H = 0, p(x,t) = 
p(x), &(x,t) = &{x), v(x,t) = 0, and S(x,t) = S(t). 

The equation of continuity (41) implies that S = —Et 


where E is a constant. Then, the Hamilton-Jacobi equa¬ 


tion (42) and the Einstein equation (44) reduce to 
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These two coupled equations determine the equilib¬ 
rium state of a spatially inhomogeneous relativistic self- 
gravitating SF/BEC. If we neglect the terms in 4>/c 2 (see 
the simplified model of Appendix |D|) , they reduce to 
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Taking the Laplacian of Eq. (B3) and substituting the 


for the density 

h 2 


result in Eq. (B4), we obtain a single differential equation 

ti 2 (Vp) 2 


h 2 


87rGm 2 \ y/p J Gm 3 


A p = p + 


8 m 2 c 2 p 


+ ■ 


2tt a s h 2 2 

3 2~P 


E 


mc- 


1 + 2 ^) P - 


In the nonrelativistic limit c —> +00, it reduces to 

47ra s fi 2 
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(B5) 


(B6) 


This returns the differential equation describing the equi¬ 
librium state of a spatially inhomogeneous nonrelativistic 
self-gravitating BEC studied in mm- 


Appendix C: The nonrelativistic limit 

In the nonrelativistic limit c —»• +00, the evolution of 
the wave function ip(x, t) in an expanding universe is de¬ 
termined by the GPP equations [lOOj : 
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These equations are equivalent to the quantum Euler- 
Poisson equations m- 
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where the pressure p is given by Eq. (45). The Eu¬ 
ler equation ( |C5| ) includes a quantum potential Q = 
—H 2 / (2ma 3 )A(^ r p) / ^fp arising from the Heisenberg un¬ 
certainty principle. In the classical limit H —> 0 (or in the 
TF limit), we recover the classical Euler-Poisson equa¬ 
tions for a barotropic fluid m- 

For a spatially homogen eous SF/BEC, the hydrody¬ 
namic equations (C3)-(C6) reduce to 
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From Eq. (C7), we find that p b oc a 


-3 


(C7) 

(C8) 

(C9) 

so the homoge¬ 


neous SF/BEC behaves as CDM. Then, Eq. (C9) implies 
a oc t 2 / 3 and p b = 1/(67rGt 2 ), which is the EdS solution. 
Combining these results with Eq. (C8), we find that the 


energy of the homogeneous SF decreases as 

9 n 

= (C10) 

Accordingly, its phase decreases as S b (t) = 
2a s h 2 /ZGm 2 t + G. To our knowledge, these re¬ 
sults have not been given before. The energy density 
and the pressure of the homogeneous SF are given by 

2 2na s ti 2 2 

e b = p b c , P b = -5— p b . (Cll) 

m a 

The case of a static Universe is recovered by taking 
a = 1 and H = 0 in the foregoing equations. For a static 
homogeneous SF, we find that S b (t) = —Et with 
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The energy of the SF is given by a sort of Einstein equa¬ 
tion “E = me 2 ” where the speed of light is replaced by 
the speed of sound. 
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Appendix D: Generalized Klein-Gordon-Poisson 
equations 

In this Appendix, we consider a simplified model in 
which we introduce the gravitational potential t) in 
the ordinary KG equation by hand, as an external poten¬ 
tial, and assume that this potential is produced by the 
SF itself via a generalized Poisson equation in which the 
source is the energy density e. This leads to the general¬ 
ized KGP equations. We then show that these equations 
can be rigorously justified from the KGE equations in 
the limit <f>/c 2 —» 0. Finally, we argue that this simpli¬ 
fied model is not sufficient to study the evolution of the 
perturbations in the linear relativistic regime. 


1. The FLRW metric and the phenomenological 
interaction Lagrangian 


We consider the FLRW metric that describes an 
isotropic and homogeneous expanding background. The 
line element in the comoving frame is 


ds 2 = g^dx^dx” = c 2 dt 2 — a{t) 2 5ijdx l dxP (Dl) 


For this metric, the d’Alembertian operator |7| writes 


Id 2 3H d _ 1 
c 2 dt 2 c 2 dt a 2 


(D2) 


In order to take the self-gravity of the SF into account, 
we introduce a Lagrangian of interaction that couples the 
gravitational potential t) to the scalar field p(x. t) 
according to 


2 

An t = (D3) 

The total Lagrangian of the system (SF + gravity) is 
given by C = + Ant- 


2. The Klein-Gordon equation 

The equation of motion resulting from the stationarity 
of the total action S = S v + Ant, obtained by writing 
SS = 0, is the KG equation 

Otyi ^ 

V + 2V(\p\ 2 ), v *+-^<S>p = 0, (D4) 

where the d’Alembertian operator is given by Eq. ( |D2[ ) 
and the gravitational potential d>(x, t) acts here as an 
external potential. For the specific SF potential ([3]), we 
obtain 


+ 



1 d 2 p 3H dp 1 
c 2 dt 2 + c 2 dt a 2 ^ 


H 2 


-v- 


S'Ka s m 

h 2 


MV = o. 


(D5) 


The energy density and the pressure, defined from the 
diagonal part of the energy-momentum tensor (131, are 
given by 


e = 


2c 2 


dip 

~dt 


+ ^lV| 2 + VM 2) , (D6) 


P = 


2 c 2 


dp 

~dt 


-^vr-vM 2 ). (D7) 


3. The generalized Poisson equation 


Eq. (D5) is the ordinary KG equation for a SF in an ex¬ 


ternal potential 4>(x, t) in an expanding background. We 
now state that $(:?, t) is actually the gravitational poten¬ 
tial produced by the SF itself. We phenomenologically 
assume that the gravitational potential is determined by 
a generalized Poisson equation of the form 


A<f> 1 

= 


(D8) 


in which the source of the gravitational potential is the 
energy density e of the SF (more precisely, its deviation 
from the homogeneous background density Cb{t))- Using 
Eq. (D6) for the energy density of a SF, and recalling the 
Friedmann equation (32), the generalized Poisson equa¬ 
tion can be written as 


A4> 

47rGa 2 


1 

2^ 


dp 


dt 


+ 


+ 


1 l2 ITT- 2 i 12 

2^^^ + 2h 2 ^ 


2na s m . 
h 2 c 2 1 


3 H 2 
8ttG' 


(D9) 


Eqs. (1D5J) and (D9) form the generalized KGP equations. 


They have been introduced in an ad hoc manner but they 
can be rigorously justified from the KGE equations (20) 
and (26) in the limit <1>/c 2 -A 0 (which, of course 
different from the nonrelativistic limit c —> +oo). 
see that t he g ravitational potential $ appears in the KG 
equation (D5) due to the cancelation of c 2 in the product 
$/c 2 x c 2 in Eq. (20). Therefore, we do not have to in¬ 


is 

We 


troduce by hand: the generalized KGP equations can 
be obtained from the KGE equations by simply neglect¬ 
ing terms of order 4>/c 2 in these equations. Similarly, the 
equations related to the generalized KGP equations (e.g. 
the generalized GPP equations, the corresponding hydro- 
dynamic equations...) can be obtained from the ones re¬ 
lated to the KGE equations by neglecting terms of order 
$/c 2 . Therefore, we do not write them explicitly (they 
are given in mb- Now, a warning is required. This 
model correctly describes the homogeneous background 
for which $ = 0 but it is not sufficient to describe the 
evolution of the perturbations in the linear relativistic 
regime because we must precisely take into account the 
terms of order $/c 2 in this regime (except, of course, in 
the nonrelativistic limit c —>■ +oo). Therefore, the use of 
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the KGE equations is mandatory to study the evolution 
of the perturbations in the relativistic regime. To show 
that the generalized KGP equations may give wrong re¬ 
sults, we derive in the next two subsections the dispersion 
relation of the perturbations in a static universe (Jeans 
problem) and the equation for the density contrast in 
an expanding universe based on the KGP equations and 
compare their expressions with those obtained from the 
KGE equations. 

Remark: We could also assume that the gravitational 
potential is determined by a Poisson equation of the form 


A<1> 

47rGa 2 ^ 


(DIO) 


in which the source of the gravitational potential is the 
pseu do re st-mass density p = \i/j\ 2 of the SF. Eqs. (D5) 
and (DIO) form the KGP equations. This approximation 
has been considered in m- However, there is an incon¬ 
sistency in coupling the relativistic KG equation (D5| to 
the classical Poisson equation (DIO). 


4. Dispersion relation in a static universe 

We consider the linear evolution of small perturbations 
in a static universe. For the generalized KGP equations, 
Eqs. (66)-(68) can be reduced to a system of two coupled 


wave equations 

1 d 2 a 
c 2 dt 2 


— A a = 1 


1 + 


E 


me* 




E \ 

d 2 8 

me 2 J 

W' 

h 2 f 

4r?i 2 

[A*- 


(Dll) 


1 d 2 8 
c 2 dt 2 


-ciAS - | 1 +f)S- 


(D12) 


Decomposing the perturbations in Fourier modes, we ob¬ 
tain the dispersion relation 


h 2 


4 m 2 c 4 

h 2 k 4 


a/ - 1 + + 


3 c 2 H 2 k 2 \ 


: 2 2ttt 2 c 2 / C 


4to 2 


+ c 2 k 2 - 4nGpb ( 1 H-y- 


= 0. (D13) 


Taking uj = 0 in Eq. (D13), we find that the Jeans length 


is determined by the equation 


+ c 2 k 2 j - 4t iGp b (1 + -£■ ) - 0. (D14) 


These equations differ from Eqs. (69) and (72) obtained 


in Sec. IV D starting from the KGE equations. 


5. Approximate equation for the density contrast 
in an expanding universe 


We consider the evolution of small perturbations in an 
expanding universe. For the generalized KGP equations, 
Eqs. (110)-(113) can be reduced to a system of two cou¬ 


pled wave equations 

d 2 S nTT dS c 2 

+ 2 H— = -4 A 6- 


H 2 


,A 2 <5+—A+, (D15) 


A4> 


4,-KGpbO 2 


= 1 + 2 -* *- 


h 2 


+ 


4m 2 c 2 a 2 
f E 


$ 

A5+-* 


E 


me 2 \ 2 me 2 


+ 1 u. 


(D16) 


Decomposing the perturbations in (spatial) Fourier 
modes, we obtain a closed equation for the density con¬ 
trast 


d 2 S k 

dt 2 


2 H 


dS k 

dt 


2 1,4 


H 2 k 


AnGpb 

i , 4-7T Gp b a 2 
1 "T k 2 c 2 


1 + 


3 c 2 


H 2 k 2 

4m 2 c 2 a 2 


~\ - ik 2 


4 = 0. 


(D17) 


In a static universe (a = 1, H = 0) , tak ing 4 = 0 
(i.e. u = 0), we can check that Eq. (D17) reduces to 


Eq. (D14) that determines the relativistic Jeans length 


derived from the generalized KGP equations. However, 
Eq. (D17) differs from Eq. (116) obtained in Sec. VIB 


starting from the KGE equations. 


Appendix E: Typical mass and scattering length of 
bosonic particles 


1. Dark matter halos 

We assume that dark matter halos are made of bosons 
in the form of BECs. We determine the mass of the 
bosons that compose dark matter halos according to 
whether they are noninteracting or self-interacting. For 
dark matter halos, we can use Newtonian gravity. The 
smallest known dark matter halo is Willman 1. It has a 
radius R = 33 pc and a mass M = 0.39 10 6 M 0 [1351 Il36l . 
We assume that this most compact halo is completely 
condensed, i.e. that it corresponds to the ground state 
(T = 0) of a self-gravitating Bose gas without radiative 
halo. 18 


18 For a value of the boson mass in the range 2.57 X 10~ 2l) eV/c 2 < 
m < 1.69 X 10 _2 eV/c 2 (see below), we have T<T C (where T c 
is the condensation temperature) for all the dark matter halos 
of the Universe (with mass M ~ 10 G — 10 11 Mq), so they can be 
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A completely condensed system of self-gravitating 
bosons (BEC) without self-interaction at T = 0 has the 
mass-radius relation MR = 9.95 h 2 /(Gm 2 ) [EH- -19 Ml- 
This gives 


m 

eV/c 2 



(El) 


Using the values of M and R corresponding to Willman 
1, we obtain a boson mass m = 2.57 1 0 - 20 eV/c 2 . 19 

A completely condensed system of self-gravitating 
bosons (BEC) with self-interaction at T = 0 in the 
TF limit has a unique radius R = n(a s h 2 /Gm 3 ) 1 ^ 2 
[571 EH Eg [67]. This gives 


("ffeH 3 (ir <-> 


Using the value of R corresponding to Willman 1, we 
obtain (fm/as) 1 / 3 (mc 2 /eV) = 0.654. In order to deter¬ 
mine the mass of the bosons, we need another relation 
(see footnote 1 in [103] '). This relation is provided by the 
constraint a/m < 1.25cm 2 /g set by the Bullet Cluster 
where a = 47ra 2 is the self-interaction cross section. 
Assuming that the bound is reached (this gives a maxi¬ 
mum bound on the mass and on the scattering length of 
the bosons) we get (a s /fm) 2 (eV/mc 2 ) = 1.7710 -8 . From 
these two constraints, we obtain m = 1.69 10 -2 eV/c 2 
and a s = 1.73 10 -5 fm. This boson mass is in agreement 
with the limit m < 1.87 eV/c 2 obtained from cosmologi¬ 
cal considerations m- 

The mass m = 2.57 1 0 20 eV/c 2 obtained for bosons 
without self-interaction gives a lower bound on the mass 
of the bosonic dark matter particle. Inversely, the 
mass to = 1.69 10 -2 eV/c 2 obtained for self-interacting 
bosons in the TF limit gives an upper bound on the 
mass of the bosonic dark matter particle. There¬ 
fore, the typical mass of the bosonic particle lies in 


considered to be at T = 0 Ez). They have a core-halo struc¬ 
ture with a solitonic core (BEC), which is a stationary solution 
of the GP equation, surrounded by a halo of scalar radiation 
in which the density decreases as r~ 3 similarly to the Navarro- 
Frank-White (NFW) |133| and Burkert ES profiles. This core¬ 
halo structure results from a process of gravitational cooling m■ 
Dwarf dark matter halos are compact objects that have just a 
solitonic core (BEC) without atmosphere. Therefore, their size 
is equal to the size of the soliton. By contrast, large dark mat¬ 
ter halos are extended objects with a core-halo structure. It is 
the radiative atmosphere that fixes the size of large dark matter 
halos. The atmosphere can be much larger than the size of the 
soliton (core). The presence of the radiative atmosphere solves 
the apparent paradox that BEC halos at T = 0 should all have 
the same radius (in the self-interacting case) or that their radius 
should decrease with their mass (in the non-interacting case), in 
contradiction with the observations [53 [58]. 

19 This result assumes that (i) Willman 1 is completely condensed 
without radiative halo, and that (ii) the observational values of 
and M ^ are accurate. More precise observational data may 
change the value of the boson mass m but its order of magnitude 
should remain the same. 


the range 2.5710~ 2 ° eV/c 2 < m < 1.6910 -2 eV/c 2 . 
The TF limit is valid for sufficiently large scattering 
lengths. An estimate of the critical scattering length 
can be obtained by substituting m = 2.57 10~ 2 ° eV/c 2 
in the relation (fm/a s ) 1 / 3 (mc 2 /eV) = 0.654. This gives 
a c = 6.0710 _59 fm. For a s < a c , the mass of the 
bosonic particle is to = 2.5710~ 2 ° eV/c 2 and, for a c < 
a s < 1.73 10~ 5 fm, the mass of the bosonic particle is 
mc 2 /eV = 0.654(a s /fm) 1 / 3 < 1.6910 -2 eV/c 2 . 


2. Boson stars 


General relativity is required to describe boson stars. 
The maximum mass of boson stars made of non¬ 
interacting bosons is M max = 0.633 hc/Gm and their 
minimum radius is R m in = 9.53 GM max /c 2 [9] . Intro¬ 
ducing scaled variables, we get 

Mmax 0 -m-ll eV / C Rmin 1 1 M max 

_ = 8 . 48 10 __ _ = i4.i_. 

(E3) 

For to = lGeV/c 2 , we obtain M max = 8.4810 20 M 0 
and Rmin = 1.20 10 -18 km. For m = 10~ lo eV/c 2 , we 
obtain M max = 0.848 M 0 and R m in = 12.0 km. For 
to = 10 -17 eV/c 2 , we obtain M max = 8.48 x 10 6 M 0 and 
Rmin = 1-20 x 10 8 km. For to = 2.57 x 10 _2 °eV/c 2 , we 
obtain M max = 3.3OxlO 9 M 0 and R rnin = 1.51xl0“ 3 pc. 
Some applications of these numerical results are given in 
the Introduction. 

The maximum mass of boson stars made of self¬ 
interacting bosons is M ma x = 0.307 hc 2 y/os/(Gm) 3 ^ 2 
and their minimum radius is R m in = 6.25 GM max /c 2 
[521 . Introducing scaled variables, we get 



We note that these results do not depend on the specific 
mass to and scattering length a s of the bosons, but only 
on the ratio m 3 /a s . For to = lGeV/c 2 and a s = 1 fm, 
we obtain M max = 1.12 M 0 and R m i n = 10.4 km. For 
(fm/a s ) 1 / 3 (?nc 2 /eV) = 0.654, we obtain M max = 6.70 x 
10 13 Mq and R m in = 20.2 pc. 


Appendix F: Numerical applications 


The complete equation governing the evolution of the 
density contrast in the matter era, Eq. (130), can be 
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written in dimensionless form as 


d 2 S 3 dS 3 
da 2 + 2 a da ^ 2a 2 


1 

i + 4- 


i + 


K 

a 

2a 
a 3 


1 + 


3(7 


<5 = 0, 


(FI) 


where we have defined k = Iz/kq, e = 
kq/kj = ( 167 T p b a 3 a 2 h 2 /Gm 4 ) 1/A , p = k h /kq = 

{^Gp b a 3 h 2 /rrfc 4 ) 1 ! 4 and v = kq/kc = 
(i kGK 2 p b a 3 /m 2 d 1 ) 1 / 4 . W e note that p = \/2v as a 
consequence of Eq. ( |131[ ) . 

It is instructive to evaluate the different quantities ap¬ 
pearing in Eq. (FI) by specifying the mass m and the 


scattering length a s of the bosons. To that purpose, we 
rewrite them in the convenient form 


= 2.16 x 10~ 10 


m 


1/2 


m 


^ e V/ c 2 


(F2) 


K '\ = 5.83 x 10’ 
m _i 


-18 


G^f er«-» 


kc 


= 1.01 X 10' 


eV/c 2 


= 6.49 x 1CT 27 , 


a = 6.19 x 10“ 


-19 ( eV/(? \ 3 / Os_ 


m ) Vfm/ 


e = 3.70 x 10 


7 ^ eV/c 2 ^ ^ ^s_y /2 


V 


= V2v = 3.00 x 10 


-17 


eV/c 2 


(F4) 


(F5) 


(F6) 


(F7) 


(F8) 


To obtain these results we have used the relation p b a 3 = 
^dm,oPo, where Qdm ,o = 0.268 is the present fraction of 
dark matter and po = 3 Hq/8ttG = 8.42 x 10 -27 kgm -3 
is the present density of the Universe calculated from 
the Hubble constant Hq = 67.11 km s _1 Mpc^ 1 obtained 
from the 2014 Planck survey. 

The numerical values of these parameters are given in 
Tables I and II for typical values of m and a s appropriate 
to dark matter halos and boson stars (see Appendix [E| . 
In the matter era (a > ai = 10 -4 ), we have 

(F9) 


Indeed, Ac/Aq = kq/kccl 3 ^ a = v/a 3 / 4 < 10 3 , 
Ac/Aj = kj/kc = vIe < 10 ~ 3 , Xq/Xh = kh/kqcl 3/4 = 
p/a 3 / 4 < 10 -3 , Xj/Xh = KH/nja 3 ^ 2 = pe/a 3 / 2 < 10 — 3 , 
and Xc/Xh = kh /Kc fl3 ^ 2 = pv/a 3 / 2 < 10 -6 . The Ta¬ 
bles make clear that relativistic effects are very weak in 
the matter era (except for perturbations with very large 
wavelengths) in agreement with the discussion given in 
the paper. 


Appendix G: The smallness of the speed of sound in 
the matter era 


The speed of sound is defined by c 2 = P'(e)c 2 , where 
e is the energy density. In the matter era, e ~ pc 2 , where 
p is the rest-mass density, so that c 2 = P'(p). For a self¬ 
interacting SF, the speed of sound is given by Eq. (59). 


It may be considered as obvious that c s <C c during the 
matter era. Actually, there is a nice way to prove this 
statement. Using p b a 3 = Qdm,oPo, and Eq. (122), we 
obtain 


c 2 __ 3a s h 2 n DM ,oHi 


2Gm 3 a 3 c 2 


(Gl) 


This quantity depends on the mass m and scattering 
length a s of the bosons only through the ratio a s /m 3 . 
It turns out that this ratio also determines the size 
of dwarf dark matter halos through the relation (see 
[571 EH [Ml [67] and Appendix |1) : 


R = 7T 


(a a h 2 \ 
\Gm 3 J 


1/2 


(G2) 


As a result, Eq. (Gl) can be rewritten as 
/3IW,o\ 1/2 H 0 R 


C-s 

C 


\ 27T 2 


,3/2. 


(G3) 


Since a > = 10 -4 in the matter era, and since the 

size of the DM halos is obviously much smaller than the 
present horizon (R c/Hq ), we conclude that c s <C c 
during the matter era. We have reached this conclusion 
without having to specify the mass and scattering length 
of the bosons. We can explicitly check on Table II that 
the condition a/a 3 <C 1, equivalent to c s c, is indeed 
satisfied for the typical values of the mass m and scatter¬ 
ing length a s of the bosons considered in the literature. 
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Physical system 

kq (m 4 ) 

Kj(m 4 ) 

£ 

o 

T 

i 

Kh{ m 4 ) 

DM halos 





a 3 = 0, m = 2.57 x 10" 2O eV/c 2 

3.46 x 10“ 2 ° 

+00 

2.60 x Hr 13 

6.49 x 10 -27 

a s = 1.73 x 10 _5 fm, m = 1.69 x 10~ 2 eV/c 2 

2.81 x 10 -11 

3.08 x 10~ 18 

1.71 x 10 s 

6.49 x 10" 27 

a s = 6.07 x 10~ 59 fm, m = 2.57 x 10~ 2O eV/c 2 

3.46 x 10“ 2 ° 

3.08 x 10~ 18 

2.60 x HT 13 

6.49 x 10" 27 

Boson stars 





a s =0, m=lx 10 9 eV/c 2 

6.83 x 10~ 6 

+00 

1.01 x 10 16 

6.49 x 10 -27 

a s ~ 0, m = 1 x 10 _1 °eV/c 2 

2.16 x 10 -15 

+00 

1.01 x 10 -3 

6.49 x 10" 27 

a s = 0, m = 1 x 10 _17 eV/c 2 

6.83 x 10 -19 

+00 

1.01 x 10" 10 

6.49 x 10" 27 

a s = lfm, m = lx 10 9 eV/c 2 

6.83 x 10~ 6 

1.84 x 10” 4 

1.01 x 10 16 

6.49 x 10" 27 


TABLE I: Values of the parameters that appear in the equation for the density contrast for different astrophysical systems. 


Physical system 

e 

a 

V 

V 

DM halos 





a s = 0, m = 2.57 x 10~ 2 °eV/c 2 

0 

0 

1.87 x HT 7 

1.33 x 10 -7 

a s = 1.73 x 10“ 5 fm, m = 1.69 x 10 _2 eV/c 2 

9.11 x 10 6 

2.22 x 10 -18 

2.31 x 10 -16 

1.65 x 10- 16 

a s = 6.07 x 10 _59 fm, m = 2.57 x 10~ 2 °eV/c 2 

1.12 x 10~ 2 

2.21 x 10 -18 

1.87 x HT 7 

1.33 x 10 -7 

Boson stars 





a s = 0, m = 1 x 10 9 eV/c 2 

0 

0 

9.49 x 10“ 22 

6.77 x 10~ 22 

a s = 0, m = 1 x 10 _1 °eV/c 2 

0 

0 

3.00 x 10" 12 

2.14 x 10- 12 

a s = 0, m = 1 x 10 _1 'eV/c 2 

0 

0 

9.49 x 10~ 9 

6.77 x 10" 9 

a s = lfm, m = lx 10 9 eV/c 2 

3.70 x 10~ 2 

6.19 x 10" 46 

9.49 x 10“ 22 

6.77 x 10~ 22 


TABLE II: Values of the parameters that appear in the equation for the density contrast for different astrophysical systems. 
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